AD-  751  785 


A  STUDY  OF  ION  TRAJECTORIES  IN  POTENTIAL 
FLOWS 


Michael  R.  Deluca,  et  al 


Ohio  State  University 


r 


Prepared  for: 


Army  Research  Office  (Durham) 


July  1972 


DISTRIBUTED  BY: 


National  Technical  Information  Service 
U.  S.  DEPARTMENT  OF  COMMERCE 

5285  Port  Royal  Road,  Springfield  Va.  22151 


BEST  AVAILABLE  CG 


UNCLASSIFIED 

Wcurtty  Cl  A  ml 


'•lion 


DOCUMENT  CONTROL  DATA  •  R  &  0 

ftwwllf  tluttlllratlon  ol  till*,  body  .1  rttlotl  mn4  In4tiinj  mtnolmll un  mini  >«  i«w«<  »h»n  lh.  own II  WWW  It  elf.ttl.dj 


I  originating  activity  (Cifpotmtm  muthot) 

The  Ohio  State  University  Research  Foundation 
1314  Kinnear  Rd. 

Columbus,  Ohio  43212 


1  REPORT  TITLE 


IS.  REPORT  SECURITY  CLASSIFICATION 

Unclassified 

>5.  QROUP  ~ 

N/A 


A  STODY  OF  ION  TRAJECTORIES  IN  POTENTIAL  FLOWS 


4  OEICRlPTivE  NOTES  (Type  ot  rmport  mnd  Inclusive  rfsfcf) 

Interim _ Technical _ 

7  AuTHORlsi  JWffti  Mswi  mlddtm  inltiml.  fair  name)  ’ 

Michael  R.  Deluca  H.  R.  Velkoff 


rm  jjj 


?A.  TOTAL  NO  OF  PAGES 

114 


7b.  no.  or  Atri 
1 


M.  CONTRACT  OR  GRANT  NO 

DA-31-124-ARO-D-246 

6.  MOJICT  NO 

20010501B700 

9rn.  OMClNATOR't  REPORT  NUMRERIH 

Technical  Report  ^11 

c. 

1D12140A14? 

d. 

9b.  other  REPORT  n 0 1 1)  (Any  o(h«r  numbtlf  thmt  mmy  bm  mmalgnmd 
thlm  rmport) 

Approved  for  public  release;  distribution 
unlimited. 

II.  luPPLCMCNTAR  V  NOTES 

\2  SPONSORING  MILITARY  AC  TIVITV 

U.S.  Army  Research  Office  -  Durham 

Box  CM,  Duke  Station 

Durham,  North  Carolina  27706 

An  analytical  investigation  of  the  motion  of  ions  in  a  fluid  was  conducted  to 
provide  insight  in  the  use  of  ions  as  a  fluid  flow  diagnostic  tool.  Assuming  small 
space  charge  density,  potential  fields  and  flows  were  considered.  The  specific 
cases  studied  were  fine  corona  wire  to  a  circular  cylinder.  Both  steady  and  non¬ 
steady  flows  were  considered,  and  the  trajectories  of  the  ions  leaving  the  corona 
wire  were  mapped. 


Ele  ctrof luldmechanlc  8 
Ele  c  tr  ohydr  odynaml  c  « 
Ele ctrogae dynamics 
Flow  Measurements 
Electrostatics 
Ions 

Fluid  Mechanics 


Unclassified 


Security  Clattificttlon 


A  STUDY  OF  ION  TRAJECTORIES  IN  POTENTIAL  FLOWS 


Michael  R.  Deluca 
H.  R.  Velkoff 

Interim  Technical  Report  #11 
Contract  DA-31-124-ARO-D-246 


U.  S.  Army  Research  office  -  Durham 
Box  CM,  Duke  Station 
Durham,  North  Carolina  27706  ‘ 

July  1972 

The  Ohio  State  University  Research  Foundation 
Columbus,  Ohio  43212 

Approved  for  public  release;  distribution 
unlimited. 


FOREWORD 


The  work  reported  herein  was  sponsored  in  part  by 
the  United  States  Aimy  Research  Office,  Durham,  under 
Contract  No.  DA-31-124-ARO-D-246.  Hie  study  presented 
herein  was  conducted  by  Mr.  Michael  R.  Deluca  in  ful¬ 
filling  the  requirements  for  a  thesis  in  his  Master  of 
Science  program  at  the  Ohio  State  University. 
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ABSTRACT 


An  analytical  investigation  of  the  motion  of  ions  in  a 
fluid  was  conducted  to  provide  insight  in  the  use  of  ions  as 
a  fluid  flow  diagnostic  tool.  Assuming  small  space  charge 
density,  potential  fields  and  flows  were  considered.  The 
specific  cases  studied  were  fine  corona  wire  to  a  circular 
cylinder.  Both  steady  and  non-steady  flows  were  considered 
and  the  trajectories  of  the  ions  leaving  the  cornona  wire 
were  mapped. 
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Definition 

"Pole",  function  of  (rx ,r2,S) 
Unit  vector  in  t] -direction 
Lift  coefficient 
Diameter  of  cylinder 
Distance  from  origin  to  wire 
Distance  t:rjm  origin  to  cylinder 


Electric  field  strength 

x -component  of  E* 

-> 

y- component  of  E 
Lift  force 
Functions 

*  >e>  >“  y 

k',£'  Dimensionless  functions 

K  Ion  mobility 

Kq  Mobility  at  original  conditions 

KpV 

M  Dimensionless  parameter  )  r,u 

P  Pressure 

P0  Pressure  at  original  conditions 

r  Polar  coordinate 

rx  Corona  wire  radius 

r2  Circular  cylinder  radius 

S  Center  spacing  between  wire  and  cylinder 

(£) 

S<p  Strouhal  Number  v  U  > 
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Ex 
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Free  stream  velocity 
Corona  wire  potential 

Drift  velocity  of  ions 

Fluid  velocity 

Ion  velocity 

Ion  velocity  in  x-direction 

Ion  velocity  in  y- direct ion 

Complex  potential  for  fluid 

Cartesian  coordinates 

Dimensionless  lengths  (x'  =  — ,  ....) 

r-> 

Angles  ions  leave  corona  wire 
Angle  ions  intersect  cylinder 
Permittivity 
flicylindrical  coordinate 
Value  of  t]  at  corona  wire 
Value  of  t]  at  cylinder 
Polar  coordinate 
Bicyllndrical  coordinate 
Circulation 
Vorticity 
Fluid  mass  density 
Electric  charge  density 
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I.  INTRODUCTION 


Duriig  the  summer  of  1968,  Dr.  Henry  Velkoff  of  the  Ohio  State 
University  performed  experiments  measuring  the  ion  flow  from  a  corona 
wire  to  a  circular  cylinder  immersed  in  a  fluid  stream  at  varying 
Reynolds  numbers.  The  experimental  data  shows  oscillations  In  the  ion 
flow  (at  approximately  the  Strouhal  frequency)  at  certain  locations  on 
the  cylinder  and  non-periodic  variation  of  the  current  at  other  loca¬ 
tions.  The  intent  of  the  tests  was  to  determine  the  efficacy  of  the 
"ion  flow"  technique  for  diagnosis  of  separation,  transition,  etc.  The 
Intent  of  this  thesis  is  to  analyze  the  same  phenomenon  from  the  stand¬ 
point  of  the  ion  trajectories. 


II.  STATEMENT  OP  PROBLEM 


Consider  a  circular  cylinder  of  radius  (r2)  in  a  fluid  stream  of 
velocity  (u) .  Directly  upstream  of  the  cylinder  is  a  corona  wire  of 
radius  (r, ).  center  spacing  (S),  and  potential  (V).  The  ions  possess  a 
mobility  (K)  and  are  leaving  the  corona  wire  at  an  angle  (a)  and  inter¬ 
secting  the  cylinder  at  an  angle  (7).  The  problem  is  complicated  by 
the  unsteady  behavior  of  the  flow  and  the  (possibly)  varying  mobility. 


This  work  will  compute  the  ion  trajectories,  taking  into  consideration 
the  unsteady  nature  of  the  flow  and  the  variation  in  the  mobility. 


III.  ANALYSIS 


Corona  Phenomena 

Corona  is  the  expression  used  to  describe  the  class  of  luminous 
phenomena  associated  with  the  current  jump  at  a  highly  stressed  elec¬ 
trode  preceding  a  spark  breakdown  of  the  gap.  Associated  with  corona 
is  an  electrode  of  small  radius.  If  such  an  electrode  is  placed  near 
one  of  low  curvature,  the  electric  field  induced  is  extremely  high  at 
the  highly  curved  electrode.  As  the  potential  difference  between  the 
electrodes  is  raised,  ionization  of  the  air  in  the  immediate  vicinity 
of  the  highly  curved  electrode  occurs;  however,  no  spark  crosses  the 
gap  because  the  field  is  too  low  farther  from  the  curved  electrode. 

The  resulting  corona  discharge  will  be  of  positive  or  negative  character 
depending  upon  the  electrode  polarity.  If  the  highly  curved  electrode 
or  wire  has  a  negative  potential,  the  positive  ions  formed  near  the 
wire  acquire  relatively  hif£i  energy  from  the  field  and  essentially  bom¬ 
bard  the  cathode  wire.  These  positive  ions,  while  forming  a  space 
charge  near  the  cathode  which  gives  the  negative  corona  its  oscillating 
characteristic,  will  produce  the  necessary  electrons  to  sustain  the  dis¬ 
charge.  The  electrons  that  are  formed  will  travel  relatively  slowly  to 
the  anode  because  of  the  reduction  in  field  strength  away  from  the  wire. 
These  electrons  will  generally  not  have  enough  energy  to  cause  further 
ionization  and  may  be  expected  to  attach  themselves  to  neutral  molecules 
of  electro -negative  gases  forming  large,  even  slower  moving,  ions.  A 
wire  of  positive  potential  will  attract  the  ionized  electrons  with  great 
intensity  causing  electron  avalanches  toward  the  wire.  This  frees  a 
positive  ion  space  charge  which  moves  toward  the  cathode  sustaining  the 
discharge  although  reducing  the  field  at  the  wire.  The  space  charge 
present  in  the  region  about  a  corona  discharge  causes  a  marked  distor¬ 
tion  of  the  electrostatic  field  (as  Poisson's  equation  predicts).  The 
fact  that  the  highly  ionized  region  constituting  the  corona  envelope 
occupies  only  a  very  small  volume  near  the  wire  permits  Ihe  theoretical 
determination  of  the  electric  field  frcan  Laplace's  equation  rather  than 
from  Poisson's  equation.  (References  2,  8). 


Electric  Field  Equations 

Appendix  B  consists  of  an  analysis  of  the  electric  field  between 
two  infinitely  long  conducting  cylinders  of  different  radii  (rx ,  r2), 
center  spacing  (s)  and  electric  potential  difference  (V).  The  results 
of  this  analysis  (which  assumed  zero  charge  density)  are 
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-2ay 

p  p  p 
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%  -  "  sinh"1^^ 
rjg  *  sinh"1  ^ 

a  =  ~  [s<-2S  (riW)  +  (rg2-ri2)*] 


Also, 

da.  *  (a2  1  r^  )'* 
dg  ■  (a2+  rg2)'4 


The  geometrical  significance  of  t]  and  a  is  shown  in  Appendix  E. 

The  corona  wire  will  be  ignored  for  the  determination  of  the  fluid 
velocity  field.  This  is  justifiable  for  two  reasons:  (l)  the  neces¬ 
sarily  small  diameter  of  the  corona  wire,  (2)  the  high  electric  fields 
in  the  Immediate  vicinity  of  the  wire  dominates  the  movement  of  the 
ions  in  that  region. 

The  diameter  of  the  corona  wire  in  most  experiments  ranges  from 
.0015"  to  .004". 
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Flow  Around  a  Circular  Cylinder 

The  motion  around  a  long  circular  cylinder  immersed  in  a  fluid 
stream  is  interesting  for  the  variety  of  changes  which  occur  with  an 
increase  in  the  Reynolds  number.  At  a  low  Reynolds  number  the  effects 
of  viscosity  are  sensible  at  large  distances  from  the  cylinder,  in 
particular  the  fluid  at  the  back  is  retarded.  At  higher  Reynolds  num¬ 
bers  two  synmetrical  standing  vortices  are  formed  at  the  back.  With 
increasing  Reynolds  numbers  these  vortices  stretch  farther  and  farther 
downstream  from  the  cylinder.  Eventually  the  standing  vortices  are 
drawn  out  to  a  considerable  length,  become  distorted,  and  break  down. 
Then  develops  the  characteristic  state  of  flow  in  which  vortices  are 
shed  alternately  find  at  regular  intervals  from  the  sides  of  the  cylinder 
with  vortex  trails  behind:  this  type  oi  flow  persists  over  a  large 
range  of  Reynolds  numbers.  The  asymmetrical  arrangement  of  the  vortices 
alters  the  pressure  distribution  around  the  cylinder.  The  eddying 
motion  has  a  definite  frequency  for  each  Reynolds  number.  Downstream 
the  vortices  assume  what  appears  to  be  a  regular  pattern.  The  vorticeB 
arrange  themselves  in  a  double  row,  in  which  each  vortex  is  opposite 
the  mid-point  of  the  interval  between  two  vortices  in  the  opposite  row. 
In  suitable  circumstances  the  trail  of  vortices  persists  for  a  consider¬ 
able  distance  downstream.  The  vortices  actually  do  not  arrange  them¬ 
selves  exactly  on  two  parallel  rows  with  a  definite  spacing  ratio.  For 
theoretical  purposes  the  system  was  considered  by  Von  Karman  to  be  com¬ 
posed  of  isolated  vortices  on  two  parallel  rows.  At  higher  Reynolds 
numbers  the  vortices  diffuse  so  rapidly  after  their  formation  that  it 
is  no  longer  possible  to  speak  of  the  formation  of  a  double  row.  At 
the  back  of  the  cylinder,  however,  vortices  continue  to  be  shed  with 
regularity  until  approximately  Reynolds  number  =  1.3  x  105.  From 
Reynolds  number  =  1.3  x  10s  to  3-5  x  ICP  there  is  loss  of  the  dominant 
periodicity  (regularity)  and  there  is  a  wide  spectrum  of  frequencies  to 
contend  with.  Above  3.5  x  106  recovery  of  the  pronounced  periodicity 
occurs.  The  range  of  Reynolds  number  1.3  x  105  to  3.5  x  10^  is  called 
the  "Critical"  and  "Post-Critical"  regime.  Above  3.5  x  10P  the  regime 
is  termed  "Transcritical. "  The  classical  description  of  the  critical 
Reynolds  number  is  (roughly):  "At  subcritical  Reynolds  numbers  the 
separation  is  laminar,  and  occurs  early,  on  the  front  of  the  cylinder. 
With  increasing  Reynolds  number,  transition  to  turbulence  in  the  bound¬ 
ary  layer  moves  ahead  of  the  laminar  separation  point,  the  now  turbulent 
boundary  layer  can  withstand  a  greater  pressure  rise,  and  separation 
moves  to  the  rear  of  the  cylinder."  Roshko  postulates  a  new  classifica¬ 
tion:  "at  subcritical  Reynolds  numbers  the  separation  is  laminar,  in 
the  supercritical  range  there  is  a  laminar  separation  bubble  followed 
by  turbulent  separation,  and  in  the  transcritical  range  the  separation 
is  purely  turbulent.  (References  4,  6,  10,  12) 

It  is  well  known  that  there  acts  on  the  cylinder  an  oscillatory 
lift  force  of  the  same  order  of  magnitude  of  the  mean  drag  and  an  oscil¬ 
latory  drag  force  superposed  on  the  mean  value.  The  oscillatory  forces 
possess  the  same  frequency  as  the  shedding  vortices  (approximately) . 
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Practically  all  the  unsteady  measurements  in  the  broad  subcritical  range 
exhibit  more  or  less  randomly  modulated  signals  at  the  Strouhal  frequency 
whether  they  represent  local  values  (e.g.  local  pressures  on  the  cylinder) 
or  integrated  values  (e.g.  forces  on  short  sections  of  the  cylinder.) 
Figure  2  shows  the  results  of  experiments  determining  the  root  mean 
square  oscillatory  lift  coefficient  as  a  function  of  the  Reynolds  number. 
(References  4,  5,  10). 

In  the  cases  of  the  motion  of  water  and  air  the  Reynolds  numbers  are 
very  large  because  of  the  very  low  viscosities  of  these  fluids.  It  is 
reasonable  to  expect  good  agreement  between  experiment  and  a  theory  in 
which  the  influence  of  viscosity  is  neglected,  i.e.  with  potential  flow 
th’fory.  The  pressure  distribution  according  to  inviscid  flow  theory 
around  a  circular  cylinder  differs  considerably  from  experiment  depending 
on  whether  the  Reynolds  number  lies  in  the  range  of  subcritical  or  super¬ 
critical  —  better  agreement  in  the  supercritical  range.  Experimental 
and  calculated  values  show  a  measure  of  agreement  on  the  front  side,  but 
at  the  rear  of  the  cylinder  the  differences  between  experiment  and  theory 
are  very  large,  and  explain  the  large  drag  force  experienced  by  a  circular 
cylinder.  The  inviscid  flow  theory  predicts  zero  drag.  Although,  gener¬ 
ally  speaking,  the  theory  of  inviscid  fluids  does  not  give  good  results 

for  drag  calculations,  the  lift  can  be  calculated  from  it  successfully. 
(Reference  13). 

It  can  be  shown  that  if  irrotational  flow  exists  within  some  portion 

of  fluid,  then  the  circulation  (k)  ,  which  is  the  line  integral  of  the 

velocity  vector  taken  around  a  closed  curve  within  a  fluid  region,  about 
any  closed  curve  is  zero  and  remains  zero;  and  the  permanence  of  irrota¬ 
tional  flow  is  established.  In  the  first  instance,  flow  is  irrotational. 
Due  to  the  lack  of  slip  at  the  boundaries,  rotation  starts  here.  Heat 
is  imparted  to  a  fluid  from  a  body  in  the  flow  in  exactly  the  same  manner 
that  vortices  diffuse  into  the  fluid.  For  very  slow  motion,  i.e.  low 
Reynolds  numbers,  heat  flows  out  in  all  directions  from  the  boundary, 
making  the  flow  rotational.  For  high  Reynolds  numbers,  the  only  fluid 
heated  would  be  in  the  narrow  layer  of  fluid  surrounding  the  body  and  in 
the  wake.  Similarly,  rotational  flow  is  confined  to  the  narrow  layer 
adjacent  to  the  boundary  and  to  the  wake.  Therefore,  the  flow  can  be 
analyzed  by  considering  the  generation  of  unsteady  vorticity  in  the  prox¬ 
imity  of  the  cylinder  and  the  resulting  feedback  on  the  velocity  and 
pressure  fields  near  the  cylinder.  (References  10,  15). 

Consider  a  cylinder  shedding  a  Von  Karman-like  vortex  trail  (Figure 
3) .  There  is  a  time  dependent  amount  of  vorticity  contained  within  the 
Curve  C  since  vorticity  in  the  wake  is  passing  through  the  curve  and  the 
rate  of  generation  of  vorticity  is  also  a  function  of  time.  For  sub¬ 
critical  and  "transcritical"  flow  the  vortex  generation  is  accurately 
periodic.  Consequently  the  amount  of  vorticity  (?)  enclosed  in  the  curve 
can  be  represented  mathematically  by 

P  =  P0  sin  wt 
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Figure  3  -  Cylinder  Shedding  Vortices 


where  £0  is  the  maximum  vorticity  of  a  single  vortex  and  u  is  the  shed¬ 
ding  frequency.  The  amount  of  circulation  along  the  curve  is  propor¬ 
tional  .to  the  amount  of  vorticity,  hence  the  circulation  ( k )  is 

k  =  k0  sin  ut 

If  we  require  the  circulation  about  the  curve  to  be  zero,  ve  must 
add  -k  to  k.  This  may  be  done  by  adding  an  isolated  vortex  at  the 
center  of  the  cylinder  with  a  vorticity  opposite  to  the  vorticity  gener¬ 
ated  at  the  cylinder  wall.  Superposing  this  isolated  vortex  to  poten¬ 
tial  flow  around  the  cylinder  enables  us  to  determine  the  feedback  effect 
of  the  unsteady  vorticity  generated  on  the  flow  upstream  of  the  cylinder. 


Fluid  Velocity  Analysis 

The  complex  potential  Wp  =  -Uz  is  for  uniform  flow  with  velocity  U 
in  the  positive  x-direction.  The  complex  potential  Wp  =  -U!r2/(z-zo)  is 
for  a  doublet  at  point  zo  with  axis  in  the  x-direction.  The  superposi¬ 
tion  of  the  uniform  flow  upon  the  doublet  yields  steady  flow  around  a 
circular  cylinder.  The  complex  potential  Wp  -  -(ijc/2jl)0n(z-zo)  is  for 
circulation  k  about  the  point  zo  in  the  negative  (clockwise  direction). 
It  can  be  shown  that  superposing  this  potential  on  the  two  previous 
potentials  still  represents  flow  around  &  circular  cylinder.  The  total 
complex  potential  is  then  (Reference  15) 

WF  =  -Uz  -  -Hif-  -  ^  Pm  (z-zo) 
f  z-zo  2TI  v  ' 

To  correspond  with  the  notation  used  in  deriving  the  electric  field 
strength  components,  we  will  call  the  radius  of  the  cylinder  (r2)  and 
the  center  (zo  =  dg).  Therefore 


wr*  (z-da> 


WF  =  -U  |"x+  iy  f  — (x-dg+iy) 
L  x-dg+  iy  J  2n 


Separating  both  sides  of  the  equation  into  real  and  imaginary  parts 
gives 


=  -U  fx^2  — K  jLtan-^-^ 
L  (x-de)^r’J  2TT  x-55 


The  velocity  components  are  computed  by  taking  minus  the  gradient  of  the 
potential  and  are 


dcpF  f  2  y2‘(x-dg)2  1  K  r 
=  "^r=UL  2  [(x-do)2+y2]2J  +  2II  |_ 


(x-^^+y2. 


V-.  .  -  M  ,  .2Ur,2  T  -yj^-1  .  JL  T _ £^2 _ 1 

dy  '  L(x-da)2+y2J  2n  L^x"d2^2+y2J 


It  will  be  found  useful  to  have  the  velocity  at  the  cylinder  sur¬ 
face  for  later  calculations  of  the  lift  force  on  the  cylinder.  Expres¬ 
sing  z  in  polar  coordinates,  the  complex  potential  is  (for  cylinder  at 
origin) 


(re*0  +  ^ 


which  may  be  separated  to  yield 


r  /  2IT 


fe  (rei0) 


=  -U^r+~-^cos  0  +  ^ 


In  polar  coordinates  the  velocity  components  are  Vr  and  Vq. 
Obviously  Vr  =  0  at  the  cylinder.  The  total  velocity  consists  of  V0. 


-u(l  +  £kf\-  JLi 

r  c^p  \  r2  /  2n  r 
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At  the  cylinder  r  =  r2.  Hence, 

vFo)r  =  r3  =  “2U  BinG 
|VF|  r.rj?  =  2U  sinO 


K 

2n r2 


K 

2JJt2 


Ion  Trajectory  Analysis 

The  ions  resulting  from  the  high  fields  of  the  corona  wire  are 
transported  from  the  wire  by  interaction  with  the  neutral  molecules  of 

the  fluid  which  flows  with  a  velocity  Vp(x,y,t).  Uncharged  molecules 

or  ions  with  zero  mobility  will  flow  at  a  velocity  of  Vp.  The  ions  of 
one  sign  coming  from  the  wire  have  associated  with  them  a  space  charge 
field  with  an  electric  charge  density  pe(y,x)  which  modifies  the  elec¬ 
tric  field  as: 


divE  (x,y) =  pe(x,y)/e 


Away  from  the  wire  we  can  assume  (because  of  the  extremely  low 
charge  density)  that  the  space  charge  density  is  negligible  and 

divE  (x,y)  =  0 

The  coupling  between  the  fluid  and  the  charges  under  the  influence 

of  the  field  gives  rise  to  a  drift  velocity  VD(x,y)  of  the  charges 
relative  to  the  fluid. 


VD(x,y)  »  KE(x,y) 

where  K  is  the  mobility  of  the  ions.  This  expression  (which  is  essen¬ 
tially  the  definition  of  mobility)  comes  from  experimental  results  which 
have  shown  that  the  velocity  of  ions  and  electrons  (in  a  dense  gas)  is 
proportional  to  the  applied  electric  field.  This  is  true  as  long  as  the 
energy  gained  by  the  ion  is  of  the  order  of  the  thermal  energy  of  gas. 

The  resulting  transport  velocity  of  the  ions  Vj(x,y,t)  is  as  follows 
VT(x,y,t)  =  Vp(x,y,t)  +  VD(x,y) 
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or 


Vj(x,y,t)  =  VF(x,y,t)  4  KE(x,y) 

For  increased  pressure  of  the  gas,  the  ion  mobility  decreases 
(over  a  large  range  of  pressure)  for  constant  temperature  approximately 
as 


where  the  subscript  (o)  denotes  the  mobility  and  pressure  at  "original 
conditions . "  Therefore 


Vr(x,y,t)  =  Vp(x,y,t)  4  K0E(x,y) 

Separating  the  velocity  equation  into  the  components  in  the  x  and 
y  directions  gives 

VIX  =  VFX  +  ■“  KoEx 


VIY  «  VFY  +  f  ^ 


The  trajectories  for  the  ion  motion  can  be  calculated  from  the 
differential  equation  that  results  from  taking  the  ratio  of  the  velocity 
components . 

dy  Po 

VIY  .  St  .  »  _  VFY  *  T 

vb  |  11  vnt*TKoEx 


The  solution  of  this  trajectory  equation  for  Bteady  and  unsteady 
flow,  constant  and  variable  mobility,  and  for  a  wide  range  of  the  system 
parameters  (as  will  be  determined  from  the  dimensionless-sizing  of  the 
differential  equation)  constitutes  the  remainder  of  this  thesis. 


x'  =  — ,  y*  =  s'  =  ,  a'  >  — ,  etc. 


Then  the  fluid  velocity  and  electric  field  components  are  written 

vra  '  u  1 1  *  +  5nJ  [(x'-vi5  *r5] 
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v  2u[  ^,(x,-dg')  * _ «_  r  .  x'-dg’  1 

^  {(x'-dg’)2  +  y,a|  2nr2  L (x'-dg1)2  +  y'2  J 

,  M  V  .  2  (  x*  (xt2+y>2-at2) _ x'-a'  j 

(’l^-1l1)r2  I  (x'2+y'‘--a'2),?  +  (2a,y1)2  (x'-a1)2  +  y'2J 


V - 3v). 


x,2+y,2+a,? 


(’>2“\)r?  I  (x'2+y,2-a,2)r  +  (2a'y*)2  (x’-a')2  +  y'2 


The  above  equations  may  be  written 

Vuy  =  U  •  f(x',y  ',dg')  +  —  •  g  (x’,y',dg») 


Vpy  =  U  •  h(x'  +  —  •  j  (x'jySdg*) 


(na'Ti  )r2 


K(x',y',a') 


%  =  7 — 15  (x,»y,»a’) 


In  general  the  pole  a  is 


a  =  “T  [s4 -2S' (r/  +  r22)  +  (r22  -  rx2)2J  * 


For  rr  <  ^  r? ,  this  pole  becomes 


a  =  ~  [s2  ■  r3Z] 

?5  L  J 


Also,  in  general,  the  distance  from  the  origin  to  the  center 
of  the  cylinder  is 

dg  =  s-dj,  =  S  -  (a2  +  r12),/4 
So,  since  a  is  of  the  order  of  V: 

de.S-a.S-i(S=-r2=) 


which  becomes 


^  =  A  [sE + tA 
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<V 


In  the  primed  coordinates  the  expressions  for  the  pole  a' 
are 


and 


a*  =  —  [S'?-l] 

2S ' 

do*  =  —  [S'^  +  l] 

2S' 

Therefore  a'  and  dg'  are  functions  solely  of  the  ratio  of  center 
distance  to  cylinder  radius.  Therefore  the  fluid  velocities  (that  are 
a  function  of  the  space  coordinates  urd  dg ' )  >  and  the  electric  field 
strength  (which  is  a  function  of  the  space  coordinates  and  a')  are  all 
functions  only  of  (besides  U,  V,  r2,  K,  (^-r^)  the  space  coordinates 
and  the  ratio  of  the  center  distance  to  the  cylinder  radius  (S')* 

So 


VpX  =  U  •  f(x',y’,S') 
VFy  =  U  •  h'(x' .y'js’) 


**  "  (n»  -ni)r2 
EV  *  (*12  -Hi)r8 


+  —  •  g’(x*  ,y'  ,s' ) 
r2 

<  “  •  J'(x’, y’,s') 
*2 

K'(x',y',s*) 

<,(x'»y,»8’) 


To  fully  understand  the  effects  of  the  nonsteady  circulation  and 
variable  mobility  on  the  ion  trajectory,  three  separate  cases  will  be 
studied:  (l)  zero  circulation,  constant  mobility;  (2)  unsteady  circula¬ 
tion,  constant  mobility;  and  (3)  zero  circulation,  variable  mobility. 
Three  separate  cases  lead  to  different  but  similar  trajectory  differen¬ 
tial  equations. 

For  case  (1)  the  trajectory  equation  becomes 


dy  dy'  U«h*  +  - — —  .  I' 

(^“ll  )*2 


dx  dx*  U  •  f  *  +  - - -  .  K' 

(n?-ni)r2 
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Dividing  the  numerator  and  denominator  by  the  free  stream  velocity 
(U)  gives 

dy'  h'  +  - — -  •  £• 

(na-ni)r2u 

KoV 

dx'  f'  +  - - - -  •  K' 


Therefore  trajectories  for  this  case  are  controlled  by  the  two 
parameters 


S', 


KoV 


(iz-ni  )r2u 

Henceforth  the  second  of  these  parameters  will  be  designated  (M) 

KoV 


M  f 


(Tl2-Tli)r2U 


For  case  (2)  the  equation  becomes 


dy' 

dx' 


h'  sin  ut  •  y  +  M  •  r 

raU _ 

f  +777  sin  ut  •  g'  +  M  •  K' 


Since  the  time  appears  in  the  trajectory  equation  case  (2)  becomes  the 
problem  of  solving  simultaneous  differential...  equations,  the  trajec¬ 
tory  equation  being  one  and  either 


or 


&.  u.r.52,tnut.g.+_ M - k> 

dt  r,  (t]2~  rji  )ra 


dy 

dt 


U  -h' 


Ko 

+  —  sin  (.ft  •  j ' 

*2 


+ 


KqV 

(T2-Tli)r2 


i' 


being  the  other  equation.  Let 


tu  tu 

2r?  ~  D 
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This  leads  to 


f.  *  sin  (—  t') 
2  dt'  raU  V  u  / 


and 


2  dt'  r,U  V  U  / 


g*  +  M  •  K* 


J'  +  M  •  /' 


As  stated  prevfously  u  is  the  shedding  frequency. 
Therefore 

^  =  2TTSt 

where  St  is  the  Strouhal  number.  So 

-  —  =  V  +  sin  (2TTSTt')  +  M  •  K' 
2  dt '  raU 


and 


-  =  h'  +  —  sin  (2rtSmt • )  •  y  +  m  • 


2  dt'  r2U 

Consequently  case  (2)  has  four  parameters 


S’ 
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&J* 


Finally  case  (3)  has  the  equation 

Po 

(jy»  h'  +  “JT  m  •  e ' 

dx'  f'  ♦  —  M  •  K* 

P 

We  must  now  determine  the  pressure  ratio.  For  an  incompressible,  steady 
flow,  the  Bernoulli  equation  is  (reference  15) 

^pVF2+P  =  |pU2  +  P0 


Therefore 


-T«  (U2-Vf2) 


2  P, 


lU 


So 


—  =  l-^U^l-f'^'2) 
Pc  2  P0 


Hence,  the  trajectory  equation  becomes 


h*  4 


M 


1  ^  ( 1  f  2  w  2\ 

dy'  _  1  2  Po  l I~f  n  > 


dx' 


f  •  +• 


M 


1-itF  («■**•*> 


.  4’ 


And  the  variable  mobility  case  is  controlled  by  the  three  parameters 


S' 


> 


M, 


1  Pl_J? 

2  P0" 


Case  (1)  and  (3)  are  ready  for  solution  now,  but  case  (2)  requires  a 
determination  of  the  proper  amount  of  circulation  for  which  to  solve  the 
equations.  As  stated  previously  inviscid  flow  theory  has  been  relatively 
successful  at  calculating  the  lift  on  bodies  inmersed  in  a  viscous  fluid 
stream.  For  steady  circulation,  the  lift  force  is  directly  proportional 
to  the  amount  of  circulation  (Streeter).  Morkovin,  Gerrard  and  Fung  all 
summarize  their  own  and  other  experiments  which  determined  the  lift  co¬ 
efficient  (oscillating)  on  a  circular  cylinder.  It  thus  seems  reason¬ 
able  to  attempt  to  find  the  proper  amount  of  circulation  from  an  inte¬ 
gration  of  the  pressure  about  the  cylinder  and  the  relation  of  this 
value  to  the  experimentally  determined  lift  coefficients.  For  nonsteady, 
incompressible  flow,  Bernoulli's  equation  is  (Reference  15) 

VF=  -  *  -  =  F  (t) 

Ot  p 

where  F(t)  is  an  arbitrary  function  of  time.  We  desire  to  find  the  lift 
force  by  integrating  at  the  cylinder,  therefore 

Vr.  =  2U  sin  n  +  — — 

F  2ur? 


and 


T’F  = 


2  Ur2cos  rp  +  —  0 

2n 


epp  =  -  2  Ur2  cos  0  t 


~  kq  sin  ut 
2n  0 
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Hence, 


Therefore, 


^pF  0  . 

1  a  *—  Kn  W  COS  (Jt 

dt  2tt  ° 


P  =  p  [F(t)  -t-  ^-2-(cos  uit)()  -  i  (2U  sin  0 


K0  4 

i  ■  —  sin 


■*)■] 


The  force  in  the  y  direction  is 


■  y  =  /  P  r,  do  sin  0 


Substituting  the  pressure  (P)  into  the  integral  gives 


Fy  =  pUk0  sin  ut  +  prsu  kq  cos  ut 


which  can  be  written 


Fy  =  p  <o  J  U2  +  (r2i*))2  [sin(u>t+7)J 


Again, 


7  =  tan”1  -y- 

7a<*> 


w  =  2n  ST  ^  =  n  ST  ^ 


N  .  v’(lM)2 

and 

jFy|  =  p,0uyrr(TTsT)- 

Which  is  interesting  because  of  the  dependence  of  the  oxciHating  lift 
force  on  the  Strouhal  number.  However  k0  is  an  unknown  quantity  and  no 
conclusions  can  be  made  regarding  the  variation  of  the  lift  force  with 
the  Strouhal  number. 

The  lift  coefficient  is  defined  by 


r  -  — * — 

L  ~  £PU‘’P 
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Substituting, 


p«0u  +  (ns<r)? 

ipU^2r, 

Hence 

Kq  =  CL 

r-jU  /l+  (TTST)2 


where  (  -!^2_  )  is  exactly  the  dimensionless  parameter  found  for  case  (2)! 
\r2U  ) 

Since  Cl  and  S-p  are  known  for  a  wide  range  of  Reynolds  numbers, 
the  equivalent  amount  of  circulation  required  can  be  calculated. 

A  calculation  of  the  drag  force  for  the  assumed  oscillating  circu¬ 
lation  showed  this  force  to  again  be  zero. 

The  corona  wire  was  neglected  in  the  fluid  velocity  analysis 
because  of  the  small  diameter  and  high  electric  field  associated  with 
the  corona  effect.  It  is  now  possible  to  determine  approximately  what 
radius  wire  is  required  for  this  assumption  to  be  valid. 

From  potential  theory,  the  maximum  fluid  velocity  about  any  circu¬ 
lar  cylinder  occurs  at  the  "shoulder"  and  is  twice  the  free  stream  veloc¬ 
ity. 


If  we  require  that  this  maximum  fluid  velocity  is  less  than  one 
percent  of  the  drift  velocity  caused  by  the  action  of  the  electric  field 
on  the  ion,  we  have  a  means  of  calculating  rr . 


Therefore 


F. 


a(Tl»-ih  ) 


(cos  hq  -  COS0g) 


and 


'I 


sinh-1  (|  .(£ 
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Hence  for  email,  i.e.  <-<  a. 


and 

Coshr^  =  ^  (e1^1  +  e  ^ )  =  JL 


So  we  can  neglect  the  cos  0  term  in  the  field  strength  equation. 
Therefore 


So 


Or 


■  ■  (Tl?-ni  )rt 

|vF|  ,  _  2U  _  2Ur,(qg-Ti1)  /rt  \ 

|Vd|\_V  5  KV  UJ 

H<  .01 - ^ -  =  .005  M 

r2  SrgCns-TjJU 


It  was  determined  that  a  suitable  minimum  (for  a  typical  experiment) 
for  M  is  1/8.  Therefore 


r,'  =  31  s  6.25  (10“4) 

r2 


In  Dr.  Velkoff *s  experiment,  the  diameter  of  the  corona  wire  was 
.0015",  the  diameter  of  the  cylinder  was  6".  This  yields  an  rx  ’  of 
2.5(10”4)  which  satisfies  the  criteria. 


TV.  DISCUSSION  OF  RESULTS 


For  case  (l),  rero  circulation  and  constant  mobility,  it  was 
decided  to  consider  three  values  of  S'  (1.5,  2,  3)>  (Dr.  Velkoff's 
experiment  used  1.5),  and  three  values  of  M  (2,  0.5,  0.125),  (Dr. 
Velkoff's  experiment  M  was  approximately  0.7).  The  ion  trajectories 
(streamlines)  for  this  case  are  represented  in  Figures  4  through  12. 

The  most  striking  feature  of  these  curves  is  the  fact  that  although  the 
trajectories  of  the  ions  are  greatly  altered  by  variation  in  the  para- 
eter  M,  the  point  at  which  the  trajectories  intersect  the  large  cylin¬ 
der  is  not  changed  very  much.  This  indicates  that  for  a  fixed  S' 
(Corona  wire  spacing),  the  average  current  reaching  the  large  cylinder 
does  not  vary  significantly  with  M.  Since  the  variation  in  M  could 
represent  a  16-fold  increase  in  the  free  stream  velocity  this  is  a 
significant  fact.  The  current  distribution  along  the  cylinder  was 
calculated  by  assuming  that  the  current  density  around  the  corona 
wire  was  uniform.  Therefore  the  current  density  at  the  cylinder  would 
be  proportioned  to  the  ratio  of  the  angle  at  which  the  ion  leaves  the 
corona  wire  to  the  angle  at  which  the  ion  "collides"  with  the  cylinder. 
So 


Where  n  denotes  the  trajectory  number,  J  is  a  dimensionless  number  pro¬ 
port  ioal  to  the  current  density,  oc  is  the  angle  at  which  the  ions  depart 
the  corona  wire,  and  y  is  the  angle  at  which  the  ions  intersect  the 
cylinder.  Figure  13  through  16  are  the  results  of  this  calculation  and 
show  clearly  the  small  effect  the  variation  in  M  haB  on  the  average 
currents  and  the  relatively  large  effect  the  variation  in  the  parameter 
S'  has  on  the  average  current.  Indeed,  there  is  only  a  14  percent 
change  in  the  maximum  current  (at  the  stagnation  point)  for  the  1600 
percent  change  in  M.  Of  course,  the  corona  phenomenon  must  be  considered. 
Since  M  is  directly  proportional  to  the  applied  voltage,  an  M  variation 
due  to  voltage  variation  would  undoubtedly  produce  a  different  number  of 
ions  at  the  corona  wire  and  a  subsequent  change  in  the  absolute  value  of 
the  currents  at  the  cylinder.  Hence  the  shape  of  the  current  of  Figure 
13  through  16  are  their  only  distinguishing  feature  for  voltage  varia¬ 
tions.  These  curves  indicate  the  currents  at  the  front  of  the  cylinder 
would  decrease  with  a  rising  free  stream  velocity.  This  contradicts 
the  experimental  evidence  from  Dr.  Velkoff's  work.  The  data  from  his 
experiment  showed  an  increase  in  the  stagnation  point  current  for  an 
increase  in  free  stream  velocity.  Dr.  Velkoff's  experiment  included 
"field  shaping"  electrodes  at  ground  potential  in  the  vicinity  of  the 
corona  wire  which  may  account  for  the  difference.  Possibly  the  higher 
drag  on  the  corona  wire  caused  a  deflection  toward  the  cylinder  which 
increased  the  current.  A  rough  calculation  showed  a  deflection  of 
0.125  inch  would  be  required  to  account  for  the  rise  in  current  measured. 
This  amount  of  deflection  is  unlikely.  Another  possible  cause  of  the 
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Figure  15.  Current  Distribution  Around  Cylinder  for 


"error"  is  the  boundary  layer  about  the  cylinder.  As  can  be  seen  from 
Figures  6  and  12,  the  trajectories  in  the  vicinity  of  the  stagnation 
point  are  "turning  up"  abruptly  near  the  cylinder.  If  the  ions  were 
"trapped"  in  the  boundary  layer  prior  to  the  time  they  "turned  up",  the 
low  velocities  in  the  boundary  layer  would  again  allow  the  electric  field 
to  dominate  the  ion  motion  and  a  smaller  (7)  would  be  obtained. 

For  case  (2),  unsteady  circulation  and  constant  mobility,  it  was 
decided  to  consider  one  value  of  S'  (2),  three  values  of  M  (2,  0.5, 
0.125),  the  maximum  value  of  the  oscillating  lift  coefficient  (C^)*  and 
a  Strouhal  number  of  0.21  (which  is  the  Strouhal  number  in  Dr.  Velkoff 's 
work.)  The  maximum  lift  coefficient  was  measured  by  Gerrard  and  is 
about  1.5  at  a  Reynolds  number  of  105.  This  Reynolds  number  corresponds 
to  the  experiment.  With  a  Strouhal  of  0.21,  the  maximum  "time"  required 
for  an  ion  to  reach  the  cylinder  is  approximately  1  of  a  cycle  of  the 
oscillating  circulation.  This  time  is  so  short  that  the  circulation  is 
virtually  constant  during  the  whole  trajectory.  Figures  17,  lB,  19  show 
the  trajectories  for  constant  positive  and  constant  negative  circulation. 
At  a  fixed  elapsed  time  the  extreme  values  of  the  ion's  angle  with 
respect  to  the  cylinder  is  shown  by  points  A  (positive  circulation)  and 
A'  (negative  circulation  on  Figure  17.  The  short  time  required  for  the 
ion  to  reach  the  cylinder  makes  the  ion  currents  sensitive  to  the 
unsteady  flow.  Consequently  it  would  be  expected  that  the  currents 
would  be  oscillating  at  the  Strouhal  frequency  and  that  these  currents 
should  be  measurable.  This  fact  was  amply  demonstrated  in  the  oscillo¬ 
scope  photographs  of  Dr.  Velkoff's  experiment.  For  purely  theoretical 
reasons,  it  was  decided  to  plot  the  ion  trajectories  for  the  case  when 
it  requires  one  cycle  ...  of  the  oscillating  circulation  for  the  ion  to 
traverse  the  gap  and  the  case  where  it  takes  "many"  cycles  for  a  traverse 
This  was  done  by  increasing  the  Strouhal  number  to  unheard  of  dimensions. 
These  trajectories  are  shown  in  Figures  20  and  21.  An  interesting  result 
was  obtained  when  an  attempt  to  find  the  time  dependent  current  at  the 
stagnation  point  was  made.  The  idea  was  to  "send"  a  large  number  of 
ions  off  the  corona  wire  at  fixed  time  intervals  which  were  much  larger 
than  the  ion  traverse  time  and  to  count  the  number  of  ions  which  reached 
an  "electrode"  at  the  cylinder  in  each  interval.  For  this  calculation, 

S’  was  2,  M  was  0.5,  CL  wt;s  1.5,  and  S<r  was  0.21.  The  cycle  was  broken 
into  twelve  intervals  and  twenty -six  "stations"  taken  around  the  corona 
wire  to  an  angle  (a)  of  30° .  This  maximum  angle  was  selected  by  consider 
ing  the  minimum  intersection  angle  (7)  which  was  produced  by  the  ion 
trajectory  for  the  appropriate  negative  circulation.  (7)  gave  the  size 
of  the  electrode  (which  was  approximately  0.4  inches  for  a  6  inch  diame¬ 
ter  cylinder) .  Ions  leaving  at  an  angle  greater  than  30°  could  not 
intersect  with  this  "electrode."  The  result  of  this  calculation  was 

Time  Interval  1  2345  6  7  8  9  10  H  12 

Ion  Count  15  11  8  7  8  11  16  20  23  24  23  20 
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Figure  20.  Ion  Trajectories  for  M  =  .125,  8*  =  2,  =  1.25 

(Ions  originating  from  common  locations  are  l80° 
out  of  phase) 


Corona 
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However,  the  above  count  Includes  only  those  Ions  reaching  the 
electrode  from  one  side  of  the  corona  wire.  Since  the  electrode  chosen 
was  symnetrical  with  regard  to  the  axis,  the  ions  leaving  the  "bottom" 
of  the  corona  wire  are  iBO  degrees  out  of  phase  with  the  "top."  Con¬ 
sequently,  to  arrive  at  the  total  amount  of  ions  reaching  the  electrode 
in  each  interval  we  must  add  the  totals  of  time  interval  (n)  with  time 
interval  (n  +  6)  (which  are  180°  out  of  phase).  If  this  is  done  all  the 
time  intervals  have  31  ions  each!  Therefore,  a  constant  current  would 
be  measured  at  the  stagnation  point  even  though  the  ion  trajectories 
are  time  varying.  Dr.  Velkoff's  data  shows  that  although  there  is  some 
variation  in  the  stagnation  point  current,  it  is  definitely  not  varying 
at  the  Strouhal  frequency  (as  so  many  currents  measured  elsewhere  are). 
An  attempt  to  calculate  the  ion  current  at  an  angle  (7)  of  54°  proved 
inconclusive.  The  "sample"  was  decreased  by  one-half  because  of  com¬ 
puter  time  limitations.  Consequently  the  total  amount  of  ions  reaching 
the  "electrode"  was  not  sufficient.  However  the  trend  toward  a  time- 
varying  current  was  seen  (as  the  experimental  data  shows) . 


For  case  (3),  zero  circulation  and  variable  mobility.  S'  was  again 

py2 

chosen  as  (0.5,  0.125)  and  1 -  took  on  two  values  (0.09,  0.33) •  The 

*=  po 

value  of  0.09  was  chosen  because  it  represents  the  maximum  for  air  for 
incompressible  flow  around  a  cylinder.  The  larger  value  (0.33)  was 
chosen  because  it  is  a  mathematical  maximum  (approximately) .  If 

py2 

\  is  larger  than  l/3,  negative  absolute  pressures  would  be  obtain¬ 
able  from  Bernoulli's  equation.  The  results  are  plotted  in  Figures  22, 

23,  and  24.  There  is  little  difference  between  Figure  22  and  Figure  12, 

m2 

of  which  the  only  difference  is  the  variable  mobility  with  **  0.09. 

Figures  23  and  24  are  interesting  but  are  probably  purely  speculative, 

PU2 

since  -  O.33  would  lead  to  compressible  flow  for  air  (and  probably 

many  other  fluids ) .  From  Figure  23  and  24  it  would  seem  probable  that 
the  problem  of  ion  trajectories  in  compressible  flow  would  be  interesting. 


Of  all  the  dimensionless  parameters  found  In  the  trajectory  analysis, 

KV 

probably  the  most  far  reaching  one  is  the  number  M  =  - - - .  As  can 

(ri2-ni  )r2U 

be  seen  from  a  comparison  of  Figures  4  through  12,  the  larger  values  of 
M  are  characterized  by  trajectories  which  are  strongly  influenced  by  the 
electric  field.  The  lower  values  of  M  have  trajectories  which  more 
nearly  simulate  the  fluid  streamlines.  M  can  be  thought  of  as  being  the 

KoV 

ratio  of  -7— - - —  to  U.  Therefore,  qualitatively,  M  is  the  approximate 

(TV’”!ll  )r2 

ratio  of  the  ion  velocity  due  to  the  electric  field  to  the  ion  velocity 
due  to  the  fluid  velocity  field. 
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V.  CONCLUSIONS 


Qualitatively  at  least,  this  analysis  verifies  the  dependence  of 
the  ion  current 3  measured  at  the  cylinder  on  the  vortex  shedding  fre¬ 
quency  (or  Strouhal  number)  characteristic  of  flow  about  a  circular 
cylinder.  The  feature  of  the  ion  trajectories  that  makes  this  so  is 
the  short  time  required  for  an  ion  to  traverse  the  corona  wire  -  cylinder 
gap  relative  to  the  period  of  the  vortex  shedding.  If  the  period  of 
the  vortex  shedding  was  small  with  respect  to  the  ion  traverse  time,  the 
trajectories  would  vary  little  and  consequently  the  ion  currents  would 
not  be  varying  sufficiently  to  be  measurable.  Figure  21  shows  this 
clearly;  especially  the  longest  trajectories.  Fortunately  Figure  21  is 
hypothetical. 

Probably  the  most  interesting  result  of  this  work  was  the  "ion- 
count"  calculation  at  the  stagnation  point.  If  anything,  it  shows  pos¬ 
sible  trouble  in  analyzing  the  ion  currents  obtained  in  experiments. 

This  calculation  shows  the  definite  periodicity  of  the  ions  reaching 
the  stagnation  point  electrode  and  also  the  fact  that  the  sum  effect  of 
the  ions  is  non-periodic .  Dr.  Velkoff's  experiment  shows  the  stagnation 
point  current  to  be  randomly  varying  (and  at  times  constant) .  It  would 
be  a  false  assumption  to  assume  the  currents  measured  at  the  stagnation 
point  are  characteristic  of  stagnation  points.  There  should  have  been 
a  similar  periodic  behavior  at  the  stagnation  point  as  there  was  at 
different  locations  in  the  front  and  rear  of  the  cylinder.  The  periodic 
current  would  have  been  manifest  at  the  stagnation  point  if  the  corona 
wire  was  rotated  to  some  angle  with  respect  to  the  wire  -  cylinder  axis. 

The  effect  of  the  variation  in  ion  mobility  for  incompressible  flow 
seems  negligible.  However,  Figures  23  and  24  indicate  that  ion  trajec¬ 
tories  for  compressible  flow,  i.e.  larger  pressure  variations,  would  be 
very  interesting. 

For  a  given  geometry  and  potential  difference  the  free  stream 
velocity  has  little  effect  on  the  average  currents  at  the  cylinder. 

This  would  indicate  that  oscillations  of  the  free  stream  velocity  would 
not  be  as  measurable  as  the  vortex  shedding  induced  oscillations.  This 
situation,  i.e.  ion  currents  measuring  one  oscillation  and  ignoring 
another  oscillation  of  the  flow  in  the  same  fluid  stream,  may  be 
beneficial  or  harmful  (it  could  lead  to  misinterpretation  of  the  current 
measurements) . 

The  average  ion  currents  are  a  relatively  sensitive  function  of  the 
spacing  of  the  corona  wire  relative  to  the  cylinder. 

These  currents  are  even  more  sensitive  than  Figure  16  indicates 
since  the  absolute  value  of  the  electric  field  (and  consequent  ioniza¬ 
tion  of  the  fluid)  would  decrease  as  the  spacing  is  increased. 
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APPENDIX  A 


FIELD  THEORY,  POISSON'S  AND  LAPLACE'S  ESJJATIONS, 
COMPLEX  POTENTIALS  AND  METRIC  COEFFICIENTS 


Preceding  page  blank 
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Consider  a  region  of  space  tj,  each  point  of  which  is  associated  one 
or  more  numbers  representing  a  physical  quantity.  The  numbers  may  speci¬ 
fy  temperature,  pressure,  velocity,  electric  field  strength,  etc.  These 
values,  for  one  kind  of  physical  quantity,  constitute  a  physical,  field. 
Therefore,  we  may  speak  of  a  thermal  field,  a  gravitational  field,  an 
electric  field.  The  mathematical  theory  of  the  subject  is  called  field 
theory,  (pp.  64-70,  Reference  9) 

— > 

A  point  P  in  euclidean  3-spa=e  is  designated  by  three  numbers,  which 
may  be  written  (Xi,X2,X3).  Suppose  there  is  associated  a  scalar  point 
function  with  each  point  in  an  arbitrary  region  tj.  This  function  may 
also  be  varying  with  time  (t). 


<P  =  cp(Xv  »X2 ,X3 ,t) 


(1) 


The  field  is  then  said  to  be  a  scalar  field.  An  example  is  the 
temperature  distribution  on  a  body. 

Or  suppose  there  is  associated  a  vector  point  function  with  each 
point  (and  time)  in  the  region  rj. 

F*=  F(Xi.  ,X2,X3,t)  (2) 


The  field  is  then  said  to  be  a  vector  field.  An  example  is  the 
electric  field  strength  in  an  electric  field. 

The  scalar  field  o.ssociated  with  a  point  is  specified  by  one  number. 
A  vector  field  associated  with  some  point  is  designated  by  three  numbers, 
Fi tF2,F3. 

F=  (Fi  »F2,F3) 

— ) 

where  Fi,F2,  and  F3  are  magnitudes  of  the  components  of  F.  Or  unit 
vectors  may  be  employed. 


F  =  aj.Fi  +  a2F2  +  aaF3 


For  any  orthogonal  coordinate  system,  the  magnitude  of  the  vector 
is 


I  F  |  =  [(Fx)2  +  (Fa)2  +  (Fa)2]'4 

Consider  any  vector  field  F  (Figure  25).  P  is  enclosed  in  a  small 


Preceding  page  blank 
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volume  At)  and  the  total  flux  of  F  through  the  surface  of  this  volume 
(at  a  fixed  time)  is 


cos  a  dA  = 


dA 


where  dA  is  a  vector  whose  magnitude  equals  the  area  dA  and  whose  direc 
tion  is  that  o£  the  _outward-drawn  normal  to  the  surface  dA.  The  angle 
a  is  between  F  and  dA. 

A  quantity  called  the  divergence  of  F*  is  defined  by  the  equation, 

$Y‘dA 

div  F  -  lim  — -  (3) 

a.,-,0  A,‘ 

This  is  a  general  definition  of  the  divergence  which  is  applicable 
to  any  coordinate  system. 

Equation  (3)  is  another  way  of  writing  the  Gauss*  theorem. 


/ 


div  F  dr)  - 


dA 


(4) 


when  S  is  the  surface  that  encloses  rj. 

Divergence  is  a  measure  of  the  strength  of  the  source  at  P*. 

Another  quantity  called  the  curl  of  F  can  be  defined  by  the  equa¬ 
tion  (again  for  a  fixed  time) 


f 


F»dS 


curl  F  a  lim 
AA->  0 


AA 


(5) 


where  the  plane  of  C  (Figure  26)  is  so  oriented  that  a  maximum  value  is 
obtained  for  the  integral,  and  where  a  is  a  unit  vector  perpendicular 
to  this  plane.  Equation  (5)  is  related  to  Stokes'  theorem  which^  asso¬ 
ciates  the  surface  integral  of  curl  F  and  the  line  integral  of  F: 


/ 


curl  F  •  dA  =  t  F  •  dS 


I 


(6) 


Fields  may  be  classified  in  terms  of  the  divergence  and  curl.  If 

div  F  -  0 
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of  all  points  in  a  region,  the  field  is  said  to  be  solenoidal  in  that 
region.  This  means  there  sue  no  sources  or  sinks  in  that  region--every 
line  of  flux  that  enters  the  region  also  leaves  the  region. 

If,  at  every  point  in  the  region, 

curl  F*  -  0 

the  field  is  said  to  be  irrotational  in  that  region. 

Consider  two  fixed  points  A  and  B  in  a  vector  field.  Figure  27. 

The  line  integral  from  A  to  B,  silong  a  path  1: 

J  F  •  dS 
(AB)i 

may  be  celled  the  potential  cpg  at  B  with  respect  to  A.  Take  another 
path  2: 


dS 


which  may  equally  well  be  said  to  define  the  potential  cpg.  Generally, 
the  two  integrals  are  not  equal,  so  that  cpg  will  not  be  a  scalar  point 
function  but  will  depend  also  on  the  path.  If  the  line  integral  depends 
on  the  path — no  scalar  potential  exists. 


e 
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For  the  case  that 


f  F  •  dS  =  f 

J(AB)i  J(A 


f  •  dS 


we  can  write 


I  F*  •  dS  -  I  f  •  dS  =  0 
-'(AB)i  J(AB)s 


Therefore 


f 

J(AB)l 


•  d S  =  -  F  •  as 


F  •  dS  +  |  F  •  dS  =  0 
J{  BA)2 


or  the  line  integral  about  a  closed  path  through  A  and  B  is 


<JF  •  dS  -  0 

And  if  this  is  true  for  any  path  in  the  field 

— > 

curl  F  =  0 

at  every  point.  Under  these  circumstances,  a  scalar  potential  cp  is 
uniquely  defined  at  any  point  B  in  the  field.  Therefore,  the  necessary 
and  sufficient  condition  for  the  existence  of  a  scalar  potential  cp  is 
that 

— > 

curl  F  =  0 

Therefore,  for  irrotational  field  at  a  fixed  time,  define 


<f'B  ‘  'PA 


f  f‘  . 


<Pr  "  9a  ‘  f  *P  -  -  I"  F*  •  dS 


dfi'  ^xT  dXl  * 


c)q>  .  dcp 

dxa  +  r—  dx3 
ox2  dx3 
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Therefore 


^2-  dx  t  +  -^2-  dx2  +  ^-2-  dxrn  -  -|Fj.dxi  +  F2CIX2  +  F3CIX3J 

And,  since  dxi,  dxp,  and  <1x3  are  independent  of  each  other, 

Fi  =  -  ,  i  -  1,  2,  3 


For  a  solenoidal  field 


Therefore, 


F  =  -  grad  q> 


div  F  =  0 


div  F  =  div  (-grad  «p)  =  -  div  grad  cp 
-  div  grad  cp  =  -  V2^  =  0 
For  the  special  case  of  Cartesian  coordinates 


v  (p  -  n  P  +  X  P  +  -s_P 


For  a  flu^.d,  the  flux  of  mass  at  a  point  (per  cross-sectional  area 
is  given  by  pVF  where  p  is  the  density  and  VF  is  the  velocity  at  some 
point.  For  no  sources  or  sinks  the  divergence  of  this  quantity  is  zero. 


div  pVF  -  0 


For  an  incompressible  fluid 


p  div  V*  *  0 


div  VF  =  0 


This  may  be  written  in  Cartesian  coordinates  as 


dU  .  dV  ,  cW  n 

5?"  ° 


where  U,  V,  W  are  the  velocity  components  in  the  x,  y,  and  z  directions, 
respectively.  This  is  the  continuity  equation  for  incompressible  fluids. 


w 


For  irrotational  fluids  the  velocity  can  be  expressed  as  the  nega¬ 
tive  of  the  gradient  of  c  potential 


VF  -  grad  <pF 

ciqpp 

U  -  T 
0  cp 

v  -  -  ry- 

The  irrotational  fluid  criteria  is  satisfied  for  inviscid  fluids. 
Therefore,  the  velocity  field  of  an  inviscid  fluid  can  be  determined 
from  solutions  of 


=  0 

(which  is  called  Laplace '  s  equation)  for  the  desired  boundary  conditions. 
Gauss'  law  for  electric  fields  is  (Reference  17) 

->  ^ 

D  .  dA  =  Q 

which  asserts  that  the  integral  of  the  normal  component  of  the  electric 
flux  density  over  any  closed  surface  S  is  equal  to  the  total  electric 
charge  enclosed  by  S. 

By  Gauss'  theorem 


div  D*  -  pe 

or  the  divergence  of  D*  at  any  point  P*  is  equal  to  the  charge  desnity  at 
that  point.  For  free  space 


Pe 


=  0 


and  therefore 


div  D  -  0 


For  an  isotropic  and  homogeneous  medium 
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Where  E  la  the  electric  field  strength  and  c  ia  the  capacitivity  (or 
permittivity) . 


Faraday' a  lav  for  electromagnetic  fields  la  (Reference  17) 


L 


dS 


0 

5t 


6  .  dA 


which  asserts  that  the  integral  of  the  tangential  component  of  the 
electric  field  strength  vector  around  any  closed  curve  C  is  equal  but 
opposite  in  sign  to  the  rate  of  change  of  the  magnetic  flux  passing 
through  any  surface  spanning  C.  If  B  (magnetic  flux  density)  is  sero 


and 


curl  E  =  0 


and  the  electric  field  is  irrotational  and  E  can  be  expressed  as  the 
gradient  of  some  potential  function 


div  D  -  div  c  E  =  pe 

Pe  ,  . 

div  i_  -  —  =  div(  -grad  cpp*) 

€ 

o  Pe 

This  equation  is  known  as  Poisson' s  equation  and  for  pe  =  0 


=  0 

which  again  is  Laplace's  equation. 

Laplace's  equation  can  be.  solved  in  a  t lumber  of  ways,  i.e. ,  mathe¬ 
matical,  graphical,  and  experimental.  The  mathematical  techniques  are 
(1)  separation  of  variables,  (2)  function/  of  a  complex  variable,  (3) 
Laplace  transform,  and  (4)  numerical  approximation  ("relaxation"  method). 
For  this  thesis  complex  variable  technique  is  the  predominant  way  in 
uhlch  this  equation  is  solved. 

An  infinite  number  of  solutions  to  the  two-dimensional  Laplace 
equation  are  easily  obtainable  from  functions  of  a  complex  variable. 


56 


Due  to  the  application  of  complex  variable  theory,  the  study  of  lnviscid 
fluid  flow  ha  a  been  greatly  expanded. 


0 

0 

0 

0 

[j 

fl 


When  x  and  y  in  the  complex  number  z  =  x  +  iy  (z  has  been  redefined) 
are  considered  variables,  then  z  is  said  to  be  a  complex  variable*  De¬ 
fining  W  as  another  complex  variable  such  that  (W  has  also  been  redefined) 
(pp.  74-77,  Reference  15). 

W  =  f(z)  =  f(x  +  iy) 

W  may  be  separated  into  its  real  Dart  and  its  imaginary  part,  called  cp 
and  i|r,  respectively 

W  =  cp(x,y,  +  it(x,y) 

where  <p  and  ♦  are  both  real  functions  of  x,y. 

The  function  f(z)  is  said  to  be  a  function  of  a  complex  variable  if 
(l)  within  some  region  there  is  one  and  only  value  of  f(z)  for  each  value 
of  z  and  that  value  is  finite  and  (2)  the  function  has  a  one-valued  de¬ 
rivative  at  each  point  within  the  region.  Within  this  region  the  func¬ 
tion  is  said  to  be  holomorphlc,  regular,  and  analytic. 


0 


Further  consideration  of  (2)  yields  relationships  that  must  be  ful¬ 
filled  by  a  function  if  it  is  analytic.  A  complex  derivative 


0 

0 

& 

I 

E 

E 

l 

E 


df  = 
dz 


lim 
&z  -»  0 


f(z  +  5z)  -  f(z) 
5z 


may  approach  its  limit  in  an  infinite  number  of  ways.  The  different 
paths  by  which  the  limit  may  be  approached  are  considered.  For  the  first 
path,  &z  is  allowed  to  approach  zero  in  the  x-direction;  i.e.,  let 
By  =  0  first,  then  take  the  limit  as  5x  approaches  zero.  This  gives 

lim  £(z  +  5z)  -  f(z)  =  iin  f (z  +  5x)  -  f(z)  =  df 

5y_>0  5x  5x->0  ^x 

Bx-*  0 

where  the  last  term  comes  from  the  second  term  which  is  the  definition 
of  a  partial  derivative.  For  the  second  path,  5z  is  allowed  to  approach 
zero  in  the  y-direction  by  letting  Bx  =  0  first;  thus 


lim 

8x-»  0 
By-*0 


f(z  +  &z)  -  f(z)  =  1  lim  f ( z  +  By)  -  f(z)  „  1  |f 
Bx  +  i  By  i  by_^0  By  i  oy 
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Since  the  derivative  rust  be  the  same  in  either  cue  if  f(s)  is  a 
function  of  a  complex  variable 


df  =  1  df 
5x  i  <5y 

However, 

t(z)  =  W  =  q  +  i* 

and  therefore, 

df  _  dg>  +  t  d^  df  ,  Ap  +  ^  q»l> 

dx  dx  dx  dy  dy  dy 

Substitution  yields 

dx  dx  i  \  dy  dy  / 

Equating  the  real  parts  and  imaginary  parts  in  each  side  of  the  equation 
yields 


dg>  _  d£  dcp  _  _  df 

dx  dy  dy  dx 

These  relations  are  called  the  Cauchy-Riemann  equations.  They  are  not 
only  necessary  but  sufficient  conditions  for  the  function  W(z)  to  be 
called  analytic. 

Differentiating  the  first  of  Equations  (7)  with  respect  to  x  and 
the  second  with  respect  to  y  and  adding  give 


+  &SL  =  o 
dx2  dy2 

which  is  the  Laplace  equation  in  two-dimensional  cartesian  coordinates. 
Therefore,  by  considering  «p  to  be  a  velocity  potential  or  an  electric 
potential,  the  real  part  of  any  function  of  a  complex  variable  is  a  pos¬ 
sible  flow  field  or  electric  field. 

Similarly,  differentiating  the  first  of  Equations  (7)  with  respect 
to  y  and  the  second  with  respect  to  x  and  subtracting  one  from  the  other 

df*  0 

dx2  dy2 
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showing  that  the  pure  imaginary  part  of  any  function  of  a  complex  varia¬ 
ble  may  also  be  the  velocity  or  electric  potential  for  a  field.  Usually 
cp  is  considered  the  potential  function,  the  if  is  called  the  stream  func¬ 
tion. 


The  functions  cp,  if  are  called  conjugate  functions;  i.e.,  the  real 
part  of  an  analytic  function  is  said  to  be  the  conjugate  of  the  imagin¬ 
ary  part.  The  curves  obtained  by 

cp(x,y)  =  constant 
>Hx,y)  =  constant 

form  an  orthogonal  system  in  the  xy-plane. 

In  the  solution  of  Laplace's  equation  by  separation  of  variables, 
the  first  step  is  to  transform  from  Cartesian  coordinates  to  a  coordinate 
system  that  fits  the  boundary  conditions.  The  work  is  expedited  by  hav¬ 
ing  a  general  method  that  allows  transformations  to  any  coordinate  system. 
Such  a  method  is  based  on  metric  coefficients  (Reference  9)» 

A  differential  length  is  expressed  in  general  orthogonal  coordinates 
(XltX3ix3)  as 

(dS)2  =  gn(dxi)2  +  g22(dx2)2  +  g33(ixa)2 
where  the  metric  coefficients  are 


where  x,  y,  and  z  are  Cartesian  coordinates. 

The  expressions  in  orthogonal  coordinates  for  the  gradient,  diverg¬ 
ence  and  Laplaces'  equation  are 


grad  q>  =  - - - 

(  6l  l ) 


where  g  =  gl  Lg2  26-33 


...  -  - 


iSL  , 

-» 

,  «>a 

b  cp 

1  .  . 

-» 

.  dtp 

1  1 
dxj. 

,  x'A 

(  822  ) 

5xa 

(1 

S33  ) ' 

*  5x3 

.  J_ 

!♦- 

i-r 

JS * 

F=]| 

dxa 

(_  822 

1 

,  833 

\j£ 

£sl!  + . 

i.r 

5<p 

|| 

1.622 

dxaj  5x3  l 

,833 

5x3_ 

II 
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The  electrical  field  around  a  line  charge  with  a  linear  charge 
density  q  is  found  by  means  of  Gauss'  law  for  electric  fields  (Gauss' 
electric  flux  theorem) ,  the  surface  of  integration  being  that  of  a 
circular  cylinder  of  radius  r  and  unit  length  coaxial  with  the  line 
charge.  If  the  charge  is  located  at  the  origin  of  coordinates  in  the 
xy-plane,  this  field  is  given  by  (Reference  14) 


The  corresponding  potential  may  be  secured  by  substituting  this  field 
into 


E  =  -  grad  cp 

and  integrating,  grad  cp  in  polar  coordinates  with  cp  a  function  of  r  only 
is 

grad  cp  =  £_ 

dr  r 


Therefore 


3.  I_  =  -  d£  r_ 
e  r2  dr  r 

and 

«p  =  -  J  (0*  r  -  t*  r0)  (8) 

The  complex  potential  function  corresponding  to  a  line  charge 
located  at  an  arbitrary  point  Zo  may  be  derived  by  means  of  the  Cauchy- 
Riemann  equations,  but  it  is  easily  written  merely  by  inspection  of 
Equation  (8).  In  polar  coordinates 

Z  =  re10 


And 


f/r  Z  =  2m  r  +  10 

Clearly  from  Equation  (8),  cp  is  the  real  part  of  -  q/e  fa  Z,  so 

W  =  q>  +  it  =  -  ^  P/r  7. 

€ 

is  the  desired  complex  function  for  a  line  charge  at  the  origin.  It 
follows  that  a  line  charge  at  Zo  has  the  complex  potential,  (Ref.  11). 

W  =  -  J  tn  (Z-Zo) 


Preceding  page  blank 
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The  function  for  n  line  charges  situated  at  Zi,  Z2, 


Kr  ■  .  »  w 


n 


W=  -  - 


OS  ^  (2-Zg) 


•••,  Zn  is  therefore 

(9) 


Let  us  superimpose  the  fields  due  to  charges  +c  at  x  =  a  and  -e  at 
x  =  -a. 


This  choice  of  the  charge  q  simplifies  the  coefficients.  From 
Equation  (9)  the  expression  for  the  complex  potential  becomes 


So 

Therefore 


W  =  -  Mz-a)  +  if.  (Z4a) 
(Z+a) 


W  =  ft! 


(Z-a) 


W  =  fo/iJ-aJ-ii') 
\x  -  a  +  iy/ 


The  denominator  of  the  natural  logarithm' s  argument  may  be  cleared  of 
the  complex  number  by  multiplying  the  numerator  and  denominator  by  the 
conjugate  of  the  denominator,  hence 


w  =  ft.  +  & +  j4  .  (x  - a  -  tel 

(x  -  a  +  iy)  (x  -  a  -  iy) 


Multiplying, 


w  -  el—  ~ &a  ^  yf.t.  |(.-:ggg) 

(x-a)2  +  y*  . 


It  can  be  shown  that 

Mu  iv)  =  Mua  +  v"2)  +  i  tan-1  ^ 

Therefore 


u 


W  =  *  for)2  +  i  tan1 

(x-a)2  +  y2 


-2ay 


x  -  a 


+  y2 


Now,  since 


W  =  ff  +  iy 


6U 


j--., 


(x-a)2  +  y2 


¥  =  tair1 


-2ay 
a2  +  y2 


tan  ¥  - 


x‘  -  a 


Multiplying  both  sides  of  the  equation  by  the  denominator  of  the  right 
side  gives 

(x2  +  y2  -  a2 ) tan  ¥  -  -2ay 

Dividing  both  sides  of  the  equation  by  tan  ¥  and  adding  the  negative  of 
the  right  side  to  both  sides  yields 

x-  +  +  *23L-  -  a2  =  0 

*  tan  ¥ 


x2  +  y2  +  —^7  »  &2 

tan  'I' 


Completing  the  square  and  adding  to  both  sides 


x2  +  v2  +  y  + 

*  tan  ¥  * 


/-a — \2  =  a2  +  ^ — a —  V 
\tan  ¥ /  \tan  ¥  / 


Therefore 


x‘  +  y  + 


tan  ¥ 


cot  ¥  = 


tan  ¥ 


1  +  cot  ¥  =  csc2¥ 


x2  +  (y  +  a  cot  ¥)2  =  a2  esc2  ¥ 


This  equation  is  a  family  of  circles  all  passing  through  the  points 
(a,o)  and  (-a,o)  with  radii  cf 


|  sin  ♦ 


and  with  centers  at 


y  =  -  a  cot  \|f 


Therefore  the  lines  of  constant  (electrical  streamlines  or  lines  of 
force)  are  nonconcentric  circles  passing  through  what  are  called  the 
poles  at  x  =  +  a,  y  =  0. 


It  also  can  be  shown  that  the  lines  of  constant  <p  are  also  non¬ 
concentric  circles.  The  radii  of  the  equipotential  circles  are  given 
by 


sinhep 


The  centers  are  at 


a  coth  q> 


The  lines  of  constant  <|<  and  constant  cp  form  an  orthogonal  net.  Hence, 
a  transformation  which  would  be  useful  for  solving  multi -cylinder  prob¬ 
lems  would  be  one  based  on  the  transformation. 

(z  -  a) 

If  we  substitute  a  coordinate  (called  r|)  for  cp  and  another  coordinate 
(called  0g)  for  we  will  have  a  transformation  which  will  produce  a 
solution  to  Laplace's  equation  so  that  the  potential  cp  is  only  a  function 
of  tj.  To  demonstrate  this  ('using  metric  coefficients)  it  is  first 
necessary  to  find  x  and  y  as  functions  of  t)  and  00. 


So  let 


n  *  19b  - (THo 


Luckily, 


Therefore, 


1  +  u  , 

fa  - - =  2  tanh 1  u 

1  -  u 


tj  +  i/9n  ~  2  tanh-1  -  =  2  coth'1  - 
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Solving  for  z. 


x  i  iy  a  eoth 


cosh 


a 


(^S) 


sinh  (-1^5) 


But 

cosh(u  +  iv)  =  cosh  u  cos  v  +  i  sinh  u  sin  v 
sinh(u  +  iv)  =  sinh  u  cos  v  +  i  cosh  u  sin  v 

Hence,  substituting 

cosh  ^  cos  ^2  +  i  sinh  3  sin  ^ 
x  +  iy  =  a  - - - sr - - - §- 

sinh  ^  cos  —  +  i  cosh  -3  sin  -2 
2  2  2  2 


The  denominator  can  be  cleared  of  the  complex  number  as  before, 


cosh  3  cos  ^2  +  i  sinh  2  sin  ^2.  sinh  2 


0B 


*  tuo  —  T  X  O  XiUi  V  £>1U  v-  OllUl  yr  COS  *  1.  COSll  J3  Sill  _ — 

x  +  iy  =  a - I - £ - 2 - El  *  ^  r 


sinh  2  cos  ^2  +  i  cosh  2  sin  ^2  sinh  2  cos  ^2  -  i  cosh  2  sin  ^2 


3 

2 


Multiplying, 


x  +  iy  -  a 


cosh  ^  sinh  £  cos2  |2  +  cosh  £  sinh  g  sin2  ^2. 

sinh2  2  cos2  ^2  +  cosh2  2  sin2  ^ 

2  2  2  2 


(sinh:  j 


sin  ^2  cos  ^2  -  cosh  ^  sin  ^2  cos  ^2 


1 


sinh2  cos2  ^2  +  cosh2  -3  sin2  22 
2  2  2  2 


0B 


Substituting  the  relationships 

cos2u  +  sin2u  =  1 
cosh2!!  -  sinh'Ti  =  1 


yields 


cosh  £  sinh  ^  -  i  sin  |2  Cos  ^2 


0B 


x  =  iy  =  a 


cosh2  ^  -  COB2  ^ 
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Substituting  ^he  relationships 


cosh  u  sinh  u  =  sinh  2u 
cosh2u  =  )>  (l  +  cosh  2u) 
cos  u  sin  u  =  p  sin  2a 
cos2u  =  (l  +  cos  2u) 


gives 


x  +  iy  -  a  T 


£  sinh  t]  -  i  \  sin  0b 


^(1  +  C08h  T])  -  j?(l  +  cos  0b) 


which  factors  to 


x  +  iy  =  a  ^  -  i  sin  0B 
cosh  T]  -  COS  0B 

Equating  real  and  imaginary  parts  on  both  sides  of  the  equation  gives 


x  =  a 


cosh  t]  -  cos  0b 


y  =  -a 


sin  0g 

cosh  t]  -  cos  0b 


The  minus  sign  can  be  dropped  since  sin(-u)  =  -sin(u)  and  cos(-u)  = 
cos(u)  and  we  can  let  0g  =  -0g  and  then  redefine  0g.  Therefore 


x  =  a 


cosh  t|  -  cos  0b 


y  =  a 


sin  0g 

cosh  t)  -  cos 


are  the  desired  transformation  functions.  This  transformation  is  called 
the  blcyllndrical  transformation.  The  name  implies  the  use  of  this 
transformation  for  solving  two-cylinder  problems.  The  inverse  trans¬ 
formation  functions  are 


T)  =  Sm 


a2)2  + 


(x-a)2  +  y* 
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and 


0 


B 


tan-1 


_ -JS3L _ 

x2  +  y*  -  a2 


For  bicylindrical  coordinates,  it  is  convenient  to  let  tj  range  from 
-  *  to  +  «  ,  which  can  be  done  by  assigning  positive  values  of  tj  to 
cylinders  for  which  x  is  positive,  and  negative  values  of  cylinders  for 
which  x  is  negative.  The  parameter  0g  represents  the  ctugle  between 
lines  drawn  from  the  poles  to  a  given  point.  Evidently  0g  =  x  repre¬ 
sents  a  line  along  the  x  axis  between  x  =  -a  and  x  =  a;  while  0g  =  0 
represents  the  remainder  of  the  x  axis.  The  portions  of  the  cylinders 
(0g  =  constant)  above  the  x  axis  are  designated  by  positive  values  of 
dg,  while  those  below  the  x  axis  are  designated  by  negative  values  of  0g. 
The  metric  coefficients  for  bicylindrical  coordinates  are  obtained  from 
(Appendix  Aj~ 


Let 


So 

and 


Tj  =  X!  ,  0B  =  xa 


2 


Differentiating, 


Simplifying, 


1  -  COS  Qg  cosh  T) 

& 


ST  (cosh  T)  -  COS  0B)J 


dy  cosh  Tj  COS  0B  -  1  _  dx 

(cosh  tj  -  cos  0B)2  Sn 


drj 


Therefore 


Also 


So 


Hence 

Substituting 


(* JLf  =  /<k\2 
\S0B/  \Si]j 


Sy  _  d  x 
dr{  cfcg 


(¥  -  (U 

®i|  ~  ®9b 


Ta(l  -  cos  0B  cosh  t])!2  f 

^  L(cosh  Tj  -  cos  0B)2  J  L 


T]  -  COB  Ogj 

Factoring  and  simplifying 


a  sin  0B  sinh  tj 
(cOSh  T)  -  cos  0B)' 


■I 


=  a 


2  1  -  2cos  0B  cosh  rj  +  cos20b  cosh  rj  +  sin20B  sinh2  tj 
(cosh  b  ■  cos  0B)4 


Substituting 


cos20b  =  1  -  sin^g 
sinh2Tj  =  cosh2r|  -  1 

2  COS2  0B  -  2COS  0B  C08h  T)  +  COSh2  T) 


gives 


a  =  a' 

bT) 


(cosh  TJ  -  COS  0g)4 


J 

3 

fl 

I] 

a 

D 

a 

a 

a 

a 

ii 

o 

o 

o 

a 

a 

o 
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Therefore 


gT| 


(cosh  n  -  cos  0R)g  _  _ af _ 

(cosh  t]  -  cos  eg )4  (cosh  tj  -  cos  00)‘ 


So 


=  g  = 


(cosh  T)  -  cos  03)' 


Also 


g  =  gfl  gr,  =  - 77 

1  (cosh  tj  -  cos  0B)4 


And 


dimension) 


Substituting 


<%B«  )%=  V 

(cosh 

ac 

n  -  cos  0B) 

I2 

,  the 

Laplacian 

is  written  (as  in 

Appendix  A  but 

i— 

f  _S*  <tSL\ 

+  A 

Lgn  ^XiJ 

dx2 

L  gsa  OX2  J  ) 

V^rp  = 

=  g-*u 

+  A  [A 

ldrl 

Lgi 

&B  L®9b 

^9bJ  ) 

But 


Hence 


For  our  problem 


V  <p  =  0 


Therefore,  dividing  by  g'%  we  are  back  to  Laplace's  equation  (but 
in  the  coordinates  rj,  03) 
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But  we  picked  tj  and  0g  to  force  the  potential  to  be  a  function  of  tj  only, 
So 


=  0 


is  the  appropriate  differential  equation,  the  solution  of  which  is  trivial 

qi  =  A  +  Brj 

Consider  two  long  parallel  metal  cylinders  with  radii  ri  and  Tz  and 
potentials  Vx  and  V2,  respectively.  They  cylinders  are  on  the  opposite 
sides  of  the  y-axie.  Designate  the  two  cylinders  by  T&.  and  T£.  Then  the 
boundary  conditions  are 


TJ  =  TX  ,  CP  »  Vl 

=  Tfe  »  <J>  =  V2 


Substitution  of  the  boundary  conditions  leads  to 


tii  -  n*? 

The  electric  field  strength  is 


{(Vam-Vin?)  +  (Vx-v2)tiJ 


where 


and 


So 


fT-=  -  grad  cp  -  -  ,aV/t  ^ 
(«t|) A  di) 


di  = 

dt) 


=  (IlzJi) 
W  -  T2/ 


(gT])*  = 


cosh  tj  -  cos  0g 


E*  -  -  at]  (>■■■  V^(cosh  tj-cos  9g) 

a\Tli*Tj2; 


It  will  be  found  useful  to  have  the  x  and  y  components  of  the  field 
strength  as  well 

dx 

and 


Ex  = 
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*2i=  gyl  .x\t 

^  (Rx?+y>'a^  + 


(2ay)2]  [(x 


-a)2  +  y2]  j 


So,  substituting,  and  letting  V  =  Vj.  -  V.> 


p  =  V  2 1  x(^V-a') 

(V’fe)  {[(x’+y^’-a2)  +  (2ay): 


[(x-a)J 


Si _ | 

+  y2]  J 


^  "  Tiv-ifeT  * 


+  V  +  a- 


[(x^+y^-a"5)  +  (2ay)2]  [(x-a)2  +  y2] 


where 


th  =  sinh'1  —  ;  i  *  1,2 
ri 


It  is  now  necessary  to  determine  the  pole  (a)  as  a  function  of  the 
radii  (rx,ra)  and  the  distance  between  the  centers  of  the  cylinders  (S). 
Call  the  distance  from  the  origin  to  the  cylinder  centers  di  and  da. 
Hence  (Reference  9) 


In  general 


di  =  a  coth  t(l  ,  dp  =  a  coth  T£ 


d  =  .  coth  ,  .  a  £2>£_n  -  «  t*1*  1  *  1 
1  Sinh  TJ  sinh  rj 


sinh  q  =  a/r 


Therefore 


d  =  a 


(a/r  y  + 


«  (a2+r2)’ 


So  the  spacing  S  between  axes  of  cylinders  is 


S  ■  d.  -  dx  =  (az+r£)%  +  (a2+rf)’4 
Solving  this  equation  for  a  gives 


a  =  2S  [S4  ’  2S?(r*  +r^  + 


7h 


Knowing  that 


siniriu  =  in  [u  +  (u2+l)*] 
The  values  for  ijl,  can  be  computed 


[>  +  ! 

(S^+  i\*l 

Lri 

/  J 

1  =  1,2 


The  sign  is  determined  by  the  side  (positive  or  negative)  of  the  y  axis 
the  cylinder  is  located.  This  is  an  arbitrary  choice.  Choose  cylinder 
No.  1  to  be  on  the  negative  side  and  cylinder  No.  2  to  be  on  the  posi¬ 
tive.  Then 


\  - 


’fe  = 


(14) 

(15) 


Hence,  given  any  two  cylinders  with  radii  (n,  t&)  and  center -center 
distance  (S),  the  electric  field  components  may  be  computed  for  any 
position.  (References  1,  9,  11,  l4). 
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NUMERICAL  SOLUTION  EQUATIONS  AND  COMPUTER  PROGRAM 
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All  the  streamline  differential  equations  (with  or  without  the 
simultaneous  velocity  equations)  were  solved  using  the  Runge-Kutta 
third  order  approximation.  (Reference  17). 

In  general,  the  simultaneous  equations, 


^=f(x,y,t)  ,  ||  =  g(x,y,t) 

at  intervals  Ax  •-  h,  can  be  solved  using 
kj.  =  f(xo,yojt0)h 
i\  =  g(xo,yo,t0)h 

k?  =  f(xo  +  l/3h,y0  +  l/3kj.,t0  +  l/3/i)h 

is  -  g(x0  +  l/3h,y0  +  l/3ki ,tQ  +  1/3 &)h 

^3  =  f(xo  +  2/3h,yo  +  2/3k/o,to  +  2/3/2  )h 

L  =  g(xo  +  2/3h,y0  +  2/3k..,to  +  2/3/2)h 

and  then  using  the  formulas 

Ay  =  J  (kx  +  3ka) 

At  =  ^  (/1  +  3/3) 

to  solve  foi  the  new  point. 

For  the  steady  state  cases  (l)  and  (3),  only  the  £&  equation  is 
necessary  for  solution.  If  the  various  increments  are  computed  in  the 
indicated  order,  each  involves  only  quantities  which  have  been  previously 
calculated . 

The  error  for  each  calculation  is  of  the  order  of  (Ax)4. 

Considering  the  longest  path  an  ion  would  be  expected  to  take  and 
allowing  the  maximum  error  for  this  path  to  be  approximately  .002,  the 
interval  for  which  the  equations  were  solved  (ax  or  Ay)  was  ,08. 

The  program  was  written  so  that  if  the  value  of  the  derivative  was 
less  than  one,  the  interval  was  taken  in  the  x-direction,  and  if  the 
derivative  was  greater  than  one,  the  interval  was  taken  in  the  y-direction. 
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This  was  done  to  minimize  the  total  error.  Also,  if  the  velocity  was 
in  the  positive  x  or  y-directlon  the  interval  was  positive  and  if  the 
velocity  was  in  the  negative  directions  the  interval  was  negative. 

List  of  Program  Symbols: 

A  -  r2 
B  = 

C  =  a 
DN  =  Ko 


JJ  =  ion  count 
El  =  Hi 
E2  =  tjo 
DM  =  M 
P  =  t 

W  =  a 

ADYDX  * 

dx' 

VIX  =  vIX 
VIY  =  Vjy 
DELX  =  AX 
DEI//  =  Ay 


YO  =  yc 
XO  =  Xo 
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Po  =  to 

Ri  =  kj 
R2  =  kg 
R3  =  k3 

51  =  li 

52  -  la 

Sa  =  I3 

RAD  =  distance  from  cylinder  center  to  ion 
ANG  =  y 


The  program  included  is  for  case  (2),  unsteady  circulation,  con¬ 
stant  mobility . 


I 

I 

I 
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-  SEXECUTE  PUFFT 

IPUFFT  99. 

DIMENSION  JJ(lf) 

DO  888  N  =1 » 1  c 
sle  jj(N)=u 
_  Dn= I • 3 1 

V  =  2j.* ( 1 0* **3 .  ) 

A=6.6666 
B=  7. 2  2 

_ C=2.78  _  _ _ 

El*-'.. 01 

E2=.AU6 

0*1780. 

_  IF(U)16u»15wtl60  _  _ 

160' ~DM=V*DN/< <E2-E1»*U*A) 

WRITE <6 *2) DM 
2  FORMAT r^OViKlw. 4  ) 

150  DO  1-c  J= 1 .26 
"  DO  30o  M=1 .6 

p=(,U79A/6.) ‘FLOAT (M-l ) 

- W=. 6666*3.1 4 16  +  T  3. 1416/150'.  AFLOAT!  J-l  J 

Y*.1*SIN(W) 

( - X  =  -C*  .  1  *C0S  TaI - 

6  ADYDXrABSI  VI  Y(  X.  Y»U.P)/VIX{X»Y»'J»P)) 

"  "  1FIADY0X-1. 120,20,30 

20  IF  t  VI  X  (  X  ,  Y  ,'J,P)  J  AO  .40,50 
40  DELX=-.Ofl 

GO  TO  130 
FtT  ~  D  EL  X  =  .0  8 

130  DD  =  VI Y ( X, Y ,u .P) /V  I  X ( X, Y »0  .P) 

-  Y0=y  -•  -  '  -  - 

XO=X 

- - —  p0  =  p 

R1=DLLX*DD 

•; - sr=Ans  i  orcx  avtxtx  .ttotp  n 

X=X0*DELX/3. 

- Y= Y0+  3 ! / 3  .  - 

P=PO+S 1 / 3 . 

- -  _DD  =  VIY  f X.Y.O.P) /Vi XIX.T  .GT.P) 

R2=DELX*DD 

- S7=TT3  ST  D'EL'XVV  \  X  C  X*  YTOTPY’' - 

<=X0+.6666*DELX 
Y=  Y0-*  •  6666* R2 
P=P0+.6666«S2 

"  "  “  DD  =  VI Y ( X,Y ,0 • P ) /Vi  X ( X, Y ,U»P) 

R3  =  DE  LX *■  DO 

a3  =  AbotDFLX7Vlx'<  X,Y,oTP‘)T 
X  =  XO  Of  L  X 
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Y«YO**25#(Rl+3.#R3> 
P*PO+  *25* I  SI  +  3,*S3 ) 


GO  TO  flu 

.. 

30  IF(VlY(X»YfU»P*)  )  60  »  60  »  70 

6^  ~'*5Ei.Y  =  -.u8' 


GO  TO  14C 
70  0ELY*.C8 

140  OD*VIY(X»Y»u»P)/VJX(XtY*UtP) 
YO*Y 
XGjX 

- po=p  — - 

R1*DELY/D0 


Sl*A9S(0FLY/VIY(X,Y.U»P)  ) -  ‘  ~  * 

Y*Y0+0ELY/3. 

XxXO+Rl/3.  .  "  * 

P*P0+Sl/3. 

1 

00= VI Y( X»Y.U*P)/VIX{X»Y»u»P) 

R2=0ELY/DD 

1  S2«ABb(DELY/VIY(X,Y»U»P) J 

Y*Y04.6666*DELY 

j  X*XO+«6666#R2 

P=PO+.6666*S2 

'  “00*V! YCX.YiUtP) /VlXCX.Y.UtP)  '  '  -  " 

R3*DEL Y/QD 

*33»ABS'(DFLY/V1  YfX'.Y.UTPj  )' 

YxYO^DELY 

- -777 

666 

X*X0+.25*<R1«-3.*R3 ) 

P=P0+.25»<61+3.»S3) 

rad=sc:rt<  <x-bi*(x-B)+y*Y)  - - 

IF(X-4. J777t3uO*30D 

I F  C RAD-A )  666»  666‘»'6  ' 

ANG*ATAN( Y/IB-XJ ) 

bib 

~TFrA  NG  -  «TST5F )  3^0  *  Hi  » b  5  5 

I F ( ANG-«996  J444»444»30u 

544" 

300 

JJfMlrJZ/iMI  +  l  *  “ 

CONTINUE 

100 

CONT I  NuE  *  .  - 

*RITE<6»77|  < JJJKJ) .KJ*1,1M 

77 

7 

FORMAT  ( 10X~»  I  5  I  4  ) 

STOP 

END 

FUNCTION  VIX<X,Y#UtP) 

A  =  6 • 6666  i 

6=7,22  i 

C  =  2.78 

El s-4.C  1 

E2  =  .  406  *  - 

0N*1 • 31 

V«2O.*(10.**3.)  -  - 

Er.X*X4  Y*Y-C*C 
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- F*E*E+4.*C*C«y*Y  •*  ~  * 

G* ( X-C ) * ( X-C )  +  V  *  Y 
DE0X=2.*( (X*E/F)-( (X-C)/G) ) 

£X*V*DEDX/<E2-£1 > 

TP=79.1*P 

VFX  =  U*(l..  +  (A*A)*(Y*Y-(X-n)*(X“B)  )/(  (  <X-B)MX-B)  +  Y*Y  ) 

r(tr-ai«(x-Bi+Y*Y)j)  -  - -  • 

l  +  1.5*U»A»Y*b!N<TP)/(6.2832#<  <  X-9  )  *  (  X-B  >  4- Y*Y  )  ) 

D=ON  "  . .  . . .  *  " 

VIX  =  VFX4[)*EX 

RETURN  . ~ 

END 

‘FUNCTION  VIYCX.Y.U.P) - 

A=6 • 6666 

-B=7 .22  - -  —  •  -  -  •—  - 

C=2 • 78 

E2S • AU6 

"ON  -  I  •  31  — - 


V=2U.«<10.*»3. ) 

‘  E*X*X  +  Y*Y-C«C 
F«E*E4A.*C*C«Y#Y 
G=(X-C)*(X-C)4y«Y 
H=X*X4Y*Y4C*C 

_DrDY  =  2V*Y^(  (H/FT-T1T7GT5 - 

EY*V*DEDy/ ( E2-E1 I 
"TT>=79.1*P  * 

VFY a  - 2 •  *  A «  A*u*  (  (  Y*(  X^3_>  \J  U  t  X~B ) «  (  X~6 )  4- y«Y  ) » I  (X“B) 
1*(X-B)4Y*Y)  )T 

1-1.5*U*A*SIM  TP )«<X-B)/ (6.28  32* ( <X-B)*( X-B)  +  Y*Y) ) 
~D  =  DM 

VI  Y  =  VFY4Di<  ey 
RETURN 
END 
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*  TWO  SIMPLE  POTENTIAL  FLOW  CASES 

I 

I 

I 

I 

1 

i 

£ 

i 

I 

I 

I 

I 

1 

85 

I 


The  general  method  developed  in  the  body  of  this  work  may  be  used 
for  any  corona  wire  geometry  or  flow  case  (the  fluid  flow  may  be  viscous 
or  inviscid).  To  demonstrate  the  nature  of  the  technique  two  simple 
cases  will  be  studied. 

The  first  case  has  an  infinite  yow  of  corona  wires  opposite  to  an 
infinite  conducting  plate.  'Tie  fluid  flow  is  directed  perpendicular  to 
the  plate. 


Figure  28.  Infinite  Row  of  Corona  Wires 
Opposite  Infinite  Plate 


The  electric  field  is  taken  to  be  constant  between  the  wires  and 
the  plate. 

Ey  -  -  ^  =  0 


The  complex  potential  for  this  flow  is  (Reference  1?) 

Wp  =  qjjr  +  Up  =  Az2 

Substituting  Z  =  x  +  iy  and  separating  the  equation  into  real  and  imagi¬ 
nary  parts  gives 


<Pp  =  Afx2-^),  \|fp  *  2Axy 
The  fluid  velocities  are  then 


'Fx 


-  2to 
dx 
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One  dimensionless  parameter  controls  this  case. 


Note  that  A  is  negative  for  the  flow  being  studied.  If 


is 


negligible  compared  to  unity  (the  maximum  value 
jectory  equation  yields  a  family  of  hyperbolas. 

than  unity,  the  equation  gives  x '  equal  to  a  constant  (the  ions  are 
travelling  perpendicular  in  straight  lines  to  the  plate) . 

The  second  case  has  two  concentric  conducting  cylinders  enclosing 
a  potential  vortex. 


of  y'  is  unity),  the  tra- 

•  jru  * 

If  J  i8  muc^  greater 


Figure  29.  Two  Cylinders  Enclosing  Vortex 

The  electric  field  for  this  case  is  radial  and  is  given  by  (Refer¬ 
ence  14) 

E>Er--‘r^f',9-0 

a 

The  complex  vortex  potential  is  (Reference  15) 


wF  =  iu/v  (z  -  zo) 


Zo  «  0,  hence 


WF  =  iuffi(z)  *  iu0f:(rei0) 


So 


Therefore 


cpp  +  icpp  =  iu?n(r)  -  u0 


cpp  =  U0 
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The  fluid  velocity  components  in  polar  coordinates  are 


Hence, 


Consequently 


1  drnF  dcpF 

~  *  “  >a  »  VFr  "  " 


vFrt  =  r  u>  VFr  =  0 


VT„  =  r  4r  =  V™  =  -  -  u 


VTr  =  ^r  =  =  JW  1 

Ir  dt  r  rir 

a 


dr 

dr 
dt  _ 

KV  1 

«■! r 

KV 

do  " 

66 

u 

Pjh  —  u 
a 

dt 

Let  r*  =  — ,  9  is  already  dimensionless 
b 


b  ^  =  J2L.  br. 
48 

a 


Integrating, 


dr'  _  KV 
r'  ”  „  b 

P,r,  ~  U 

a 


KV 

PvT'  =  — 9  *  C’ 

u 

a 


which  may  be  written 


r»  «  c  e^'n  u 
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which  is  obviously  a  spiral.  This  case  also  has  one  dimensionless 
parameter 


KV 
„  b 

in  -  n 

a 

Substituting  initial  conditions  of  r*  =  2:  and  0  =  0  gives 


r’ 


e 


KV 

in  -  u 
a 


9 


as  the  trajectory  equation. 
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e-  C 


Figure  30.  Eicylindrical  Coordinates.  The  lines  r)  ■  constant  are 
circles  with  axes  on  the  x-axis.  The  lines  ■  const) 
cure  portions  of  circles  with  axes  on  the  y  axis. 
(Reference  9)- 


All  the  trajectory  analyses  undertaken  in  this  work  have  been 
primarily  the  solution  of  the  differential  equation 

dy  a  VFy  +  *ev 
dx  Vp*x  +  K  Ex 

Another  method  of  determining  the  ion  trajectories  for  steady  fl'nr 
which  will  lead  to  closed  solutions  is  simply  adding  the  stream  function 
of  the  fluid  flow  (\ty)  to  the  product  of  the  mobility  times  the  stream 
function  of  the  electric  il  field  (iE) .  This  sue  yields  a  combined 

fluid-electrical  stream  function  (^)  which  may  then  be  solved  for  the 
trajectories. 

♦fE  “  +  Kl^E 

In  general,  I  would  state  that  this  method  is  valid  for  steady 
state,  constant  mobility  ion  flow  problems. 
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APPENDIX  G 


ABSTRACT  OF  EXPERIMENTAL  TESTS  OF  FLUID  FLOW  DIAGNOSTICS 
USING  ELECTROSTATIC  CHARGES 
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Studies  were  initiated  into  the  possible  use  of  charged  particles 
in  a  gas  stream  as  a  means  of  studying  the  nature  of  gas  flow.  Two 
basic  approaches  were  followed.  In  one  case,  ions  were  generated  elec¬ 
trostatically  by  means  of  a  corona  discharge  from  a  suitable  point  or 
line  source.  In  the  second  case,  the  charges  were  obtained  either  from 
the  natural  dust  or  water  particles  in  the  air,  or  dust  particles  were 
deliberately  added  to  the  air  stream.  The  charges  were  collected  at 
selected  positions  downstream  on  the  particular  aerodynamic  shape  being 
studied. 


Experimental  Study  With  Ions 

Three  configurations  have  been  studied.  They  include  a  six-inch- 
diameter  cylinder  mounted  transverse  to  the  air  stream,  a  three-foot- 
long  sharp-edged  flat  plate,  and  an  NACA  0012  airfoil  of  0.216  m  chord 
length,  all  of  which  were  mounted  in  the  U.  S.  Army  AMRDL  7  x  10  tunnel 
at  Ames  Research  Center  and  tested  at  various  tunnel  speeds.  Most  of 
the  runs  were  made  using  a  0.0038  cm  diameter  corona  wire  placed  a  few 
centimeters  upstream  of  the  surface.  Electric  field  shaping  electrodes 
adjacent  to  the  wire  were  grounded.  Strips  of  aluminum  1.23  cm  wide, 
foil  spaced  at  0.32  cm  were  mounted  on  each  surface.  An  end  view  of 
one  configuration  is  shown  in  Figure  31.  The  filter  is  used  to  elimi¬ 
nate  background  60  Hz  noise. 


Field  shaping 
y^clcctroacs 

V 

\ 


Voltage 

Figure  31.  Ion- Flow  Measurement  Scheme 
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The  corona  currents  ranged  from  2  uA  to  30  \jA  and  the  applied 
potential  used  to  get  the  corona  ranged  from  6  to  l6  kV  which  depended 
mainly  on  the  configuration  of  the  electrodes.  Essentially  the  same 
results  were  achieved  regardless  of  the  currents  used. 

Figure  32  shows  a  typical  oscilloscope  trace  using  the  cylinder 
at  a  tunnel  q  of  one.  The  upper  traces  are  taken  from  the  resistor 
connected  to  the  most  forward  position  on  the  cylinder  and  the  lower 
grace  is  from  the  most  aft  position.  The  regularity  of  the  upper  traces 
indicates  that  the  ion  currents  tend  to  follow  the  basic  oscillating 
flow  about  a  cylinder  in  cross  flow.  Such  flow  oscillations  are  the 
result  of  the  periodic  shedding  of  vortices  from  a  cylinder.  At  the 
rear  of  the  cylinder  the  flow  has  largely  separated  and  close  to  the 
cylinder  is  essentially  a  "dead-water"  region.  Examination  of  the  lower 
trace  reveals  that  the  current  oscillation  is  reduced  greatly  or  elimi¬ 
nated.  The  frequency  of  the  oscillation  is  l6  Hz  which  is  of  the  correct 
order  for  the  Strouhal  frequency  for  this  cylinder.  Thus,  the  ion  tech¬ 
nique  in  this  case  seems  to  respond  properly  to  the  characteristics  of 
the  flow  field. 


Figure  32.  Circular  Cylinder  Trace 


Next,  a  flat  plate  was  placed  directly  aft  of  the  cylinder  at  the 
centerline  and  extended  0.46  m  rearward.  Such  surfaces  act  to  eliminate 
or  reduce  the  periodic  shedding.  Figure  33  shows  the  results  using  the 
ion  technique.  The  upper  trace  is  from  the  front  of  the  cylinder  and  it 
can  be  seen  that  the  ion  technique  indeed  does  show  greatly  reduced 
oscillation. 


Figure  33.  Circular  Cylinder  With  Aft  Plate 


The  three- foot -long  sharp-edged  flat  plate  was  mounted  in  the 
tunnel  parallel  to  the  windstream.  The  electrodes  were  located  on 
the  upper  surface.  The  top  trace  of  Figure  34  illustrates  an  oscilla¬ 
tion  in  flow  0.10  m  aft  of  the  leading  edge.  These  oscillations  may 
indicate  the  presence  of  a  separation  bubble  near  the  leading  edge. 

The  next  trace  from  a  position  3.8  cm  aft  of  the  first  trace,  although 
it  shows  excursions,  does  not  show  oscillations.  The  lowest  trace 
further  back  on  the  plate  shows  a  relatively  smooth  appearance.  The 
ion  technique  seems  to  give  a  proper  indication  of  the  flow  for  this 
case. 
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Figure  3*+.  Sharp- Edged  Flat  Plate 


Tests  nerct  were  run  on  the  airfoil  using  the  corona  wire  and  field¬ 
shaping  electrodes.  Current  collecting  strips  were  located  over  the 
entire  airfoil.  At  a  tunnel  q  of  one  (8.2  m/sec)  and  an  angle  of  attack 
of  0°,  all  the  traces  from  positions  around  the  airfoil  were  smooth  and 
regular  except  at  one  location.  The  upper  trace  of  Figure  35  is  10.8  cm 
from  the  nose  of  the  airfoil.  The  next  two  traces  are  1.6  cm  and  3.2  cm 
further  back  on  the  0.216  m  chord  airfoil.  It  can  readily  be  seen  that 
a  significant  oscillation  exists  within  a  limited  region  on  the  airfoil, 
indicating  the  presence  of  a  separation  zone.  Thus,  the  ion  technique 
appears  to  be  able  to  discriminate  local  phenomenon.  All  the  foregoing 
ion-flow  tests  were  taken  with  a  filter  setting  which  filtered  out  fre¬ 
quencies  about  30  Hz.  This  is  the  reason  why  no  high  frequency  oscilla¬ 
tions  appear  in  the  data. 

Another  test  using  the  airfoil  set  at  an  angle  of  attack  of  17°  and 
a  tunnel  q  of  10  (V  =  26  m/sec)  was  .  *n  with  the  filter  set  to  pass  fre¬ 
quencies  in  the  150-2000  Hz  range.  Figure  36  illustrates  violent  oscilla¬ 
tions  of  the  ion  currents  at  two  pot>:  tions  near  the  midchord  region  of 
the  upper  surface.  This  action  is  expected  since  the  airfoil  is  in  full 
stall.  The  lower  surface  traces  show  relatively  random  behavior  with 
peak  amplitudes  of  approximately  1-5^  of  the  upper  surface  trace  ampli¬ 
tudes. 


Figure  36.  Airfoil  at  17®  Angle 


Charged  Particles  Tests 

The  next  phase  consisted  of  using  the  charged  particles  carried  in 
the  wind  tunnel  stream  as  the  source  of  charge.  The  configuration  used 
was  the  cylinder  with  the  rear  splitter  plate.  The  electrode  at  the  rear 
edge  of  the  plate  was  monitored.  Large  oscillations  can  be  noted  from 
the  trace  in  Figure  37.  The  flow  in  this  region  contains  large  scale 
turbulence  and  such  large  scale  oscillations  are  expected.  If  corona 
currents  are  used  with  this  test  configuration,  similar  oscillations  are 
seen  but  their  magnitude  is  several  times  greater. 


Figure  37.  Charged  Particle  Traces 


